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ABSTRACT 

Since the discovery of the Kuiper Belt, broadband surface colors were thoroughly studied as a first approximation to the object 
reflectivity spectra. Visible colors (B VRI) have proven to be a reasonable proxy for real spectra, which are rather linear in this range. 
In contrast, near-IR colors (JHK bands) could be misleading when absorption features of ices are present in the spectra. Although the 
physical and chemical information provided by colors are rather limited, broadband photometry remains the best tool for establishing 
the bulk surface properties of KBOs and Centaurs. In this work, we explore for the first time general, recurrent effects in the study of 
visible colors that could affect the interpretation of the scientific results: i) how a correlation could be missed or weakened as a result 
of the data error bars, ii) the “risk” of missing an existing trend because of low sampling, and the possibility of making quantified 
predictions on the sample size needed to detect a trend at a given significance level — assuming the sample is unbiased, iii) the use 
of partial correlations to distinguish the mutual effect of two or more (physical) parameters, and iv) the sensitivity of the “reddening 
line” tool to the central wavelength of the filters used. To illustrate and apply these new tools, we have compiled the visible colors 
and orbital parameters of about 370 objects available in the literature — assumed, by default, as unbiased samples — and carried 
out a traditional analysis per dynamical family. Our results show in particular how a) data error bars impose a limit on the detectable 
correlations regardless of sample size and that therefore, once that limit is achieved, it is important to diminish the error bars, but it 
is pointless to enlarge the sampling with the same or larger errors; b) almost all dynamical families still require larger samplings to 
ensure the detection of correlations stronger than ±0.5, that is, correlations that may explain ~25% or more of the color variability; 
c) the correlation strength between (V-R) vs. (R-I) is systematically lower than the one between (B-V) vs. (V-R) and is not related 
with error-bar differences between these colors; d) it is statistically equivalent to use any of the different flavors of orbital excitation or 
collisional velocity parameters regarding the famous color-inclination correlation among classical KBOs — which no longer appears 
to be a strong correlation — whereas the inclination and Tisserand parameter relative to Neptune cannot be separated from one 
another; and e) classical KBOs are the only dynamical family that shows neither (B-V) vs. (V-R) nor (V-R) vs. (R-I) correlations. It 
therefore is the family with the most unpredictable visible surface reflectivities. 

Key words. Kuiper belt: general — Methods: data analysis — Methods: statistical — Techniques: photometric 


1. Introduction 

Kuiper Belt objects (KBOs), also known as Trans-Neptunian 
objects (TNOs), are a large population of icy bodies orbiting 
the Sun with semi-major axes greater than that of Neptune. 
They are considered as the least altered objects in the solar sys¬ 
tem. Their existence had been hypothesized by Leonard| ( 1930| l, 
Edgeworth ( |1943[ |1949| l, and Kuiper ( |1951| l — and the expres¬ 
sion Edgeworth-Kuiper Belt is also used in some works — but 
they were only observed in 1992 by Jewitt & Luu ( 1993| l. Cur¬ 
rently, there are more than 1700 objects identified, ~ 4.5% are 
binary or multiple systems, and there are probably ~ 200000 
objects larger then 100 km in diameter ( |Petit et al.|2011[ ). 

KBOs subdivide into several dynamical families, although 
there is no strict definition regarding the boundaries between 
most of them (see [Gladman et al.|[2008[ for a review). Cen¬ 
taurs are, presumably, former KBOs that evolved into inner or¬ 
bits, that is, both their semi-major axes and perihelia lie be¬ 
tween the orbits of Jupiter and Neptune (|Eernandez| 19801 |Lev-| 


ison & Duncan|[TW7 1 — although when using a less strict dy¬ 
namical definition of Centaurs, some may have an Oort Cloud 


origin ( Emel’yanenko et al.|2013[ Eouchard et al.|20T4) l. Every 
~ 125 yr a new Centaur escapes from the Kuiper Belt into a 
short-lived chaotic orbit, with an average half-live of ~ 2.75 Myr 
— although some may last for millions of years — before be¬ 
coming a Jupiter-family comet (JEC) or being re-injected into 
the outer solar system (e.g., [Horner et al.|2004| l. Given this close 
link. Centaurs have usually been studied as one more family of 
KBOs. 


Since the first analysis by Luu & Jewitt| ( 1996) 1 of a sample of 
three Centaurs and nine KBOs, any relation between the surface 
properties observed at different wavelengths has been searched 
for, such as different colors. These relations might be among the 
surface properties themselves or between the surface properties 
and the orbital parameters. Although correlation analysis does 
not allow us to predict specific values of unmeasured proper¬ 
ties by knowing others, nor to draw conclusions about cause and 
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effect relationships, it is the most elementary tool for detecting 
and quantifying which variables or properties covary with each 
other, which provides a first insight into the phenomena we are 
observing. 

Many works reported trends, patterns, and correlations be¬ 
tween different surface colors and between those and the orbital 


properties of these objects (see review Doressoundiram et al. 
2008|l. While several teams had their own surveys — for ex 


Tegler & Romanishin 1998 2000 2003[ 

Tegler et al. 

2003|l, 

the ESO Large Program on Centaurs and I 

fNOs (e.g., Delsanti 

et al.|[2001 

Boehnhardt et al.[[2002 [Delsanti et al.|[2004| [Peix- 

inho et al. 

2004 Delsanti et al.[[2006[l, the Meudon Multicolor 

Survey (2MS), (e.g., [Barucci et al.|1999 2000| [Doressoundiram 

etal.|2001| 

2002| [2005b |2007|l, and the Second ESO Large Pro- 

gram (e.g.. 

DeMeo et al.|2009| Perna et al.|2010 and references 


therein) — focusing on the analysis of their own datasets for 
methodological reasons, |Hainaut & Delsaiiti] ( |2002 1 made the 
first thorough analysis of a large compilation of all published 
data, known as the MBOSS database, which has recently been 
updated ( [Hainaut et al.|201^ . 

Although the physical and chemical information provided by 
the visible colors of these objects is rather limited, the analysis 
of visible colors alone has lead to very rich debates on the Kuiper 
Belt science. 

Given the narrow temperature differences in the Kuiper Belt 
region, KBOs were initially presumed to have a primordial ho¬ 
mogeneous composition, and the first modeling works attempted 
to explain the visible surface color diversity as a result of a 
competition between a reddening space-weather effect of irra¬ 
diation on an ice-rich surface layer and a bluishing collisional 
resurfacing by buried non-irradiated layers (Luu & Jewitt 1996 1 . 
A strong correlation of color vs. orbital inclination detected 
among the dynamic family of classical KBOs ( [Trujillo & Brown| 
2002| l led to a further development of collisional resurfacing 
models, accounting for differently irradiated subsurface layers 
(|Gil-Hutton]|2002|l, or for collision-triggered cometary activity 


and size-dependent resurfacing (Delsanti et al. 2004|l. While 


these models had difficulties in explaining the detection, or 
non-detection, of several surface properties and correlations be¬ 
tween these and the orbital properties (e.g., Jewitt & Luu|200l] 


Thebault & Dores soundi ram|2003[|Delsanti et al. 2004| l, the dy 

namical model of |Gomesfp003| l suggested that high-inclination 
KBOs formed closer to the Sun migrate afterward into their cur¬ 
rent orbits. This opened a new window for the possibility of orig¬ 
inal compositional differences among KBOs (e.g., ITegler et al.l 
[ 2 ^ . 

Simultaneously, a debate on the global existence, or absence, 
of only two distinct surface color groups of Centaurs and KBOs 
( |Tegler & Romanishin||199^ |20()3| ) led to the conclusion that 
only Centaurs seemed to be separated into two distinct groups, 
whereas KBOs presented a somewhat continuous distribution of 
colors ( [Peixinho et al.|[2003| ), and current data suggest that the 
existence of two distinct groups of objects is related to their sizes 
and not to their dynamical families ( [Peixinho et al.|20T^ Fraser 
|& Brown|2012| l. 

Dynamical models indicate that the giant planets have mi¬ 
grated throughout the history of the solar system, sculpting the 
structure of the current Kuiper Belt (e.g., [Levison et al.||20lT] 
Batygin et al.||20lT and references therein). Although there is 


still a debate on the details of this process, it became more plau¬ 
sible that the Kuiper Belt is composed by a mixture of bodies 
formed at distinct heliocentric distances and, therefore, their dis¬ 
tinct surface properties would be primarily caused by primordial 
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compositional differences (e.g., [Brown et al.|2011[ ). The most re¬ 
cent analysis of albedos vs. visible colors of Centaurs and KBOs 
identifies a surface color and albedo separation that is also ev¬ 
idence for such primordial compositional differences ( [Lacerda[ 
et al. 2014|l. Furthermore, latest laboratory works show that it 


is possible to reproduce a wide range of surface colors of Cen¬ 
taurs and KBOs with an appropriate combination of initial albe¬ 
dos, collisional evolution, and space weathering, implying that 
it is unlikely that a straightforward dependance between the dy¬ 
namical properties and the surface properties exists ( [Kanucho^ 
et al.[2012|. Thus, it comes as no surprise that the two taxonomy 


schemes proposed, so far, identify several groups with distinct 
surface properties (Barucci et al.|2005| Perna et al.|2010[ [Fraser[ 


& Brown[2012[ Dalle Ore et al.[20T3 1, failing to detect any clear 

dependance between the dynamical properties and the surface 
properties, although the proportions of taxonomic families vary 
among the dynamical families. 

With the growing interest on detailed spectroscopy studies in 
the near-IR, where several ice absorption bands can be detected 
, as well as with the existence of an apparently large database of 
visible colors of our objects, it might be perceived as irrelevant to 
continue photometric studies of KBOs and their progeny. Flow- 
ever, spectroscopic studies are still limited to fewer than 100 ob¬ 
jects (i.e., ~ 5% of known objects) and are extremely biased to¬ 
ward only the brightest objects, whereas visually all detected ob¬ 
jects might be studied photometrically, and, as aforementioned, 
it is currently known that smaller (fainter) objects do possess a 
peculiar bimodality of surfaces that is not seen among the larger 
(brighter) objects. The question is, nonetheless, pertinent, as in 
theory one does not need to study all the population but only a 
representative sample of it. The problem is that the more subtle 
the effect one seeks, the larger the sample required. 

In this work we exclusively focus on studying the correla¬ 
tions between the visible surface colors and the orbital parame¬ 
ters for each traditional dynamical family of KBOs and Centaurs 
of our own compilation of visible colors available in the litera¬ 
ture, assumed to be unbiased samples of each family, by default. 
We analyze and discuss for the first time,what we might be miss¬ 
ing as a function of the available sample sizes for each dynamical 
family, and also as a function of the observational error-bars. For 
an analysis of the (statistical) distributions of the surface proper¬ 
ties of each dynamical family see[Hainaut et al.[([2012[l. 


2. Data sample 

Orbital elements were gathered from the Asteroid Orbital Ele¬ 
ments Database, astorb.dat[^ maintained by the Lowell Ob¬ 
servatory based on astrometric observations by the Minor Planet 
Center. Orbital inclinations i are relative to the ecliptic and not 
to the so-called Kuiper Belt plane. 


We used the classification scheme suggested by Lykawka & 
Mukai ( 2007[l, including their analysis of objects located in the 


mean motion resonances (MMR) with Neptune, following a ten- 
step algorithm: 


1. in 3 : 2 MMR with Neptune ^ Plutino 

2. in other MMR with Neptune other resonant 

3. q < aj ^ not analysed 

4. q > aj A a < ^ Centaur 

5. aj < q < A a ^ scattered disk object (SDO) 

6. qm < q ^ 31AU ^ scattered disk object (SDO) 

7. q>AQAUAa>A^AU^ detached KBO (DKBO) 

* ftp://ftp.lowell.edu/pub/elgb/astorb.dat.gz, with epoch 20120420 
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8. 31AU < ^ < 40At/ scattered or detached KBO (SD- 
KBO) 

9. i < 5° A{[q > 31 AU A (37At/ < a < 40At/)] y [q > 
38At/ A (42 AU < a < 48At/)] ) ^ cold classical KBO 
(cCKBO) 

10. i ^ 5° A q ^ 31 AU A (31 AU < a < 48 At/) hot classical 
KBO (hCKBO). 


We are aware that there are more complex classification 
schemes that may be more rigorous, but the boundaries between 
families do not change significantly. We chose this one for its 
computational simplicity, which provides a dynamical classifica¬ 
tion for all objects in our sample, whereas dynamical simulation- 
based classifications are only available for a limited number of 
objects (e.g., |Gladman et al.|2008| l. 


The so-called Haumea family (|Brown et al.|2007||Snodgrass 
let al.||2010|l — (136108) Haumea, (24835) lOOSSMjs, (55636) 
2 OO 2 TX 300 , (86047) I 999 OY 3 , (19308) 1996X066, (120178) 
2 OO 3 OP 32 , (145453) 2 OO 5 RR 43 , 2003SQ317, 2003UZii7, and 
2 OO 5 CB 79 — are extremely peculiar objects, therefore we re¬ 
moved them from the statistical analysis. We also removed the 
two objects in retrograde orbits ( 2 OO 8 KV 42 , 2 OO 8 YB 3 ). For a 
total of 366 objects in our database, we carried out a global anal¬ 
ysis, that is, all objects except Haumea family and retrograde 
orbits, of n = 354 objects, calling it All*; another analysis also 
removed Centaurs, which we call KBOs*; a third analysis ac¬ 
cording to dynamical family; and finally, an analysis for which 
we combined some families. 

We compiled the visible colors (BVRI broadband filters or 
equivalent filters transformed to that system) available for 366 
objects (see Table]^ and fully trusted their published error-bars. 
For each input paper and each object or observation, we com¬ 
puted the reflectance spectrum using Eq. (3) from |Delsanti et ah] 
( 20011 1, when two or more filters were available. The resulting 
spectra were manually checked and obviously deviating data 
from a given filter were removed from the dataset. Color indexes 
were computed within one given epoch, which led to colors ob¬ 
tained from simultaneous photometry (the different bands were 
observed over a maximum time-span of two hours). Then the 
average colors indexes and their one-sigma errors from different 
papers and epochs were computed for each object using Eqs. (1) 
and (2) from |Hainaut & Delsantij f2002| l, which provided more 
accurate estimates when multiple measurements are available. 
The absolute magnitude in R band (//*) was computed for each 
object and epoch whenever an R-band magnitude was available, 
using Hr = R-5 log(rA), where R is the R-band magnitude and r 
and A the helio- and geocentric distances at the time of observa¬ 
tions, respectively. Different values for a given object were also 
averaged using Eqs. (1) and (2) from |Hainaut & Delsanti| ( 2^002| l. 
We did not correct for any phase effect, as already discussed in 
[Peixinho et al.| ( |2012| l; see Table 


3. Statistical methods 

3.1. Statistical tests 

- Spearman-rank correlation accounting for data error- 
bars: We analyzed the possible existence of correlations be¬ 
tween surface colors of Centaurs and KBOs and orbital pa¬ 
rameters using the Spearman-rank correlation coefficient p 
( |Spearman|1904| l. The significance of p was estimated using 
t = p[(n-2)/(l -p^)]*^^ , which obeys approximately a Stu¬ 
dent t-distribution with n - 2 degrees of freedom, where n is 
the number of data-points ([Kendall et al.|1939]l. 


To take into account the effect of observational errors, we 
considered each data-point not as fixed, that is, exact, but 
as a Gaussian distribution centered on the measured value 
with a standard deviation equal to the corresponding error- 
bar. Then, we extracted 1000 new random data-sets in which 
each data-point followed its own Gaussian probability dis¬ 
tribution. Since the Spearman p -distribution is not Gaus¬ 
sian, we used the Eisher transformation ^(p) = ^ ln(|i^) = 
argtanh(p) — Eisher ( 1915|l —, determined the mean {(p), 
that is, the mok probable ^p value, the H-cr^Z-cr^, and trans¬ 
formed it back to p-values using p^ = tanh(^p), obtaining the 
most probable correlation value (p^), and the shortest inter¬ 
val containing 68.2% of the distribution, which is the equiv¬ 
alent of a Gaussian +lcr limits, which in this case are not 
necessarily symmetrical. We denote this as -Hcr^Z-crg (see 
[Peixinho et al.|2004 for more details). 

This method yields the confidence interval of the most proba¬ 
ble p of our sample, but does not yield the confidence interval 
for the correlation value of the parent population. Eor that we 
made 1000 bootstrap extractions from our data-points and 
computed the corresponding {pt) and +cri,l-crh, analogously 
to the previous case (e.g., [Efron & Tibshirani|19931l. 

Usually, the ±crb interval dominates over ±ae, but a more 
accurate estimation of the confidence interval for the correla¬ 
tion of the parent population is given by quadratically adding 
V(o-c + cPA, which, after reconversion gives us the best esti¬ 


mate for the 68.2% confidence interval: 

- Partial correlation 

When dealing with multivariate variables, we constantly face 
the same question: given three variables A, B, and C, is 
the correlation observed between the two variables A and 
B (Pab) due to the association of A and B with the third 
variable, C, and not to the direct association between them? 
Using the partial correlation, we may eliminate the suspected 
effects of this third variable C on the relation between A and 
B by making C somehow constant, denoted by Pab.c- This is 
achieved by 


Pab-RacPbc 

PAB.C = , (1) 

V(1 “Pac^)(1 - Pbc^) 

If Pab.c does not differ significantly from pab , then the as¬ 
sociation between A and B is most probably direct and not 
masked by their correlations with C. To take into account 
the effects of data error-bars, we have used the same pro¬ 
cess as previously described for the Spearman-rank correla¬ 
tion. The significance of pab.c is estimated analogously to 
the case of p , but with one less degree of freedom, that is, 
using t - Pab.c [(« - 3)/(l -P 4 Rr)]'^^ — see, for instance, 
[MacklinldEi^ . 


3.2. Assessing the risk of missing existing correlations 

Searching for correlations means testing a null hypothesis (Hq), 
that is, variables X and Y do not show evidence for correlation, 
against an alternative hypothesis (ffi), where variables X and Y 
show evidence for correlation. This evidence is measured by the 
significance level (S L), confidence level (CL — 1 - S L), or p- 
value (= SL). The threshold p-value below which we reject Hq is 
called a. The lower the significance, the lower the probability of 
incorrectly rejecting the null hypothesis (type I error). To dimin¬ 
ish the odds of finding false correlations, we simply adopted a 
stricter criterion for rejecting Hq, for instance, instead of accept¬ 
ing an a = 1% probability of falsely finding a correlation, we 
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accept only a - 0.1%, or even lower. The drawback is that this 
increases the chances of not detecting correlations that might in 
fact exist, meaning that we increase the risk, j3, of missing some 
real correlations (type II error), or decrease the (1 -yS) power 
to find them. If we wish to diminish this risk, we must either 
increase the a threshold with which we reject the null hypothe¬ 
sis or increase the sampling. We have analyzed the sample size 
requirements, n, to detect a correlation at some given signifi¬ 
cance level (higher or equal to a) and to only miss a fi propor¬ 
tion. Since we wish to compare the case where the correlation is 
indeed p Q with the case of p = 0, from|Modarres|([T996|l and 


Bonett & Wright (20001 we obtain 


4(zi-„/2+zi-p)"(l+pV2) 

■ 

where Zi-a/ 2 , and z\-p are the quantiles of a normal distribution 
for our threshold a/2 — since we are using two-tailed tests — 
and risk /?, respectively, from which we compute the smallest 
integer not lower than n, that is, [n], ceiling of n. For example, to 
see evidence for a correlation at a 2.5 cr level, that is, a - 0.0124 
and Z\~ai 2 - 2.5, in a sample extracted from a parent population 
which we knew to possess p = 0.6 accepting only a /I = 0.1 
probability (risk) of missing it, that is, Z\-p - 1.2816, we would 
need a sample of n > 39. If we accept a higher risk of missing 
the detection, for example, (3 - 0.2 {z\-p - 0.8416), we would 
need only n > 31. 



3.3. Degraded correlations 


When taking into account the data error-bars, the possible cor¬ 
relations between the X and Y variables will diminish or de¬ 
grade. That is, taking each data-point not as an exact value but 
as a Gaussian probability distribution centered on the measured 
value diminishes the possible association strength between the 
two variables. The larger the error-bars, the stronger this effect. 

The error-bars of surface colors of Centaurs and KBOs tend 
to become larger with increasing wavelength (e.g., from visible 
B - V to the near-infrared H - K). It might be possible that two 
different pairs of color indexes, Xi vs. Yi and X 2 vs. Y 2 , corre¬ 
late with equal strength (e.g., B - F vs. V - R and V - R vs. 
R - I) — if we precisely knew the full population, for instance. 
However, not only are we estimating the populations’ ps from 
a limited sample, but if the error-bars of (X 2 , Y 2 ) are larger than 
those of (2fi, Yi) we might also obtain a |p 2 | < |pi| or even lose 
evidence for correlation on (X 2 , T 2 ). To analyze this possibility, 
we performed the following study: 


1. We generated a parent population of 1000 (X, Y) pairs with a 
desired correlation value. Each variable follows a Gaussian 
distribution, and their values are renormalized into the inter¬ 
val [0,1] (the mean will become 0.5) by two distinct meth¬ 
ods: a) the lowest value will correspond to zero, the highest 
to one; b) the mean value minus 3 cr will correspond to zero, 
the mean value plus 3 cr to one (with this renormalization 
some points may fall outside the [0,1] interval). 

2. We extracted from that parent population a sample of n el¬ 
ements. We added Gaussian error bars to the Y elements 
and generated 1000 new samples from the n elements by 
the same process as we used for the real data-samples (see 
Sect.|3.1|l. The approximate sampling error given by [Bonett] 
& Wright| ( |2000l l, 


a2^(l+p2/2)/(n-3). 
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(3) 


was added quadratically to the dispersion of the Fisher- 
transformed simulations — instead of using bootstraps ex¬ 
tractions — before reconverting them to obtain a final . 

3. We iterated the two previous steps for p = 0.8,0.7,0.6, and 
0.5, each with samplings of n = 15,20,30,50,75, and 100, 
and varyed the observational error-bars from cTohs = (Jo - 
0.00 to 0.30 in steps of 0.02. 

4. We repeated these cycles with X and Y following a random 
uniform distribution instead of a Gaussian. 

First, the simulations show that drawing the sample from 
a Gaussian distribution or from a uniform distribution, either 
renormalized into the interval [0,1] by the methods la) or lb), 
has little to no influence on the degradation of p values and their 
error-bars, with increasing simulated observational errors cTo . 
Variations are on the order of 0.01 or lower. Hereafter, the anal¬ 
ysis refers to the simulations drawn from a Gaussian distribution 
that was renormalized using method la). 

Second, we find that the effects of the (simulated) observa¬ 
tional errors on X and on Y, which we denote by <Tox and <Toy, 
are equivalent to performing the simulations with additional ob¬ 
servational errors only on the X variable of cto - Uox + cToY, 
adding no errors to the Y variable. For example, having a series 
of both X and Y values with error-bars of 0.1 has an equivalent 
effect on the correlation values as having a series of X values 
with error-bars of 0.2 and Y values with no errors. This result is 
most convenient and simplifies the rest of our analysis. 

Third, we find, empirically, that for each initial p value its 
degradation with increasing observational errors cTo is indepen¬ 
dent of the sample size n. The estimated 68.2% error interval 
on the measured changes with n because it is largely de¬ 

pendent on it (see Eq. |^, but the diminishing of (p) shows no 
evident dependence on n. 

Using a Fevenberg-Marquardt algorithm ( jLevenbei^j 1 944[ 
|Marquardt||1963| ), we find that third-degree polynomials fitted 
with great accuracy (all cases with residuals rms < 0.0063) the 


detectable p functions denoted as p* = /(p, cTo)'. 



' 1.000 - 0.950cro - 10.544cr2 -h 25.3490-^ 

p = 1.0 


0.900 - 0.497cro - 8.668cr2 h- 18.3870-^ 

II 

0 


0.800 - 0.309cro - 7.255cr2 h- 14.4350-^ 

00 

d 

II 

P* = ■ 

0.700 - 0.176cro - 6.373cr2 h- 12.2990-^ 

p^O.l 


0.600 - 0.129cro - 4.322cr2 h- l.A6\crl 

p - 0.6 


0.500 - 0.05130-0 - 2.8170-2-H 4.3730-3 

p = 0.5 


(4) 

An immediate consequence of the p degradation with larger 
(ToS is the magnitude decrease of the detectable p*. That is, even 
if we had a parent population such that X -Y, hence p = 1, the 
fact that a sample has error-bars will impose a limit on the de¬ 
tectable correlation, regardless of the sample size n. In practice, 
from the correlation analysis point of view, this means that it is 
very important to diminish the observational errors, but rather 
pointless to enlarge the sampling with the same or larger error- 
bars if the elements in the sample are already representative of 
the whole population, of course, and the significance level of the 
detection has reached the threshold we desire. 

Another consequence is that increasing /? increases the risk 
of missing an existing correlation. Therefore, to keep p constant 
(for a theoretically detectable correlation), we need to increase 
the sampling n. We return to the example in Sect.|3.2| We know 
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the parent population has p - 0.6 and we need n > 31 to detect 
evidence for a correlation at a 2.5 cr level with a risk /? = 0.2 of 
missing it. Suppose now that each X data-point has an observa¬ 
tional error cTox-O-^ and each Y data-point also has ctoy^O.I. As 
we saw earlier, this is equivalent to transferring all the error-bars 
to the X variable, such that <To - cTox + ctoy - 0.2, and to assume 
Y with no errors. (Note that we take X and Y as normalized to 
[0,1] so this error represents 10% of the full range of values on 
X and 10% on Y.) From Eq. we see that with such observa¬ 
tional errors the detectable correlation drops to p* = 0.46, and 
now, from Eq. to detect evidence for a correlation at the same 
2.5 cr significance level, we will need « > 53 data-points, that 
is, 70% more data-points (n > 67 for (3 - 0.1). In this work, we 
are able to quantitatively predict the sample size needed for the 
cases studied below. 


3.4. Correction for the false-discovery rate of the correlations 

Since for each sample or subsample we performed about m - 
150 correlation tests, the chances of finding evidence for a cor¬ 
relation increase, even if there is no correlation at all (see |Mill^ 
1981 for a review). Eor example, suppose we reject the null hy¬ 


pothesis of no correlation using a threshold a = 0.0124 (2.5 cr 
level). There is a 1.24% probability of obtaining evidence for a 
correlation just by chance, even if it is not there. This means that 
performing 150 tests will probably lead to one or two false dis¬ 
coveries. To solve this problem, we used the false-discovery rate 
(EDR) procedure as proposed by |Benjamini & Hochberg ( 1995| l. 
That is, instead of ensuring that the probability of erroneously 
rejecting even one of the null hypothesis is always lower than a 
certain a — using the Bonferroni procedure of rejecting Hq only 
when its p-value is < ajm ( |Bonferroni|1935| l —, we ensure that 
the proportion of incorrectly rejected Hqs, or false discoveries, 
is a. 

The Benjamini-Hochberg procedure works as follows: 


1. Order the m p-values from all tests as pi < p 2 < ■■■ < Pm, 

2. find the largest integer k such that k = max{i: pi < ^a), 

3. reject all the null hypotheses whose p-value is < pk. 


This procedure ensures that EDR < a. That is, if we use a - 
0.0124 (2.5 cr) at most 1.24% of all the postive results (i.e., all the 
rejected null hypotheses) may be false positives. Eor example, if 
we simultaneously test 150 hypotheses and reject fewer than 30 
of them, or in other words, we have fewer than 30 discoveries, 
it is unlikely that we make even one false discovery (at 2.5 cr, 
naturally). 


4. Reddening line 


Delsanti et al.| ( |20()T] ) introduced the concept of drawing a line in 
the color-color plots indicating the location of objects that proba¬ 
bly possess linear reflectivity spectra. As discussed in detail later 
by Delsanti et al.| (20041, a linear reflectivity spectrum does not 
result in a linear relation between colors, and vice versa. This 
line has been widely used to infer the linearity or non-linearity of 
object spectra. However, each telescope has its own set of filters, 
and the differences between the effective central wavelength for 
the same filters among the different sets has a significant effect 
on the location of this reddening line. This effect can be strong 
for the visible BVR colors, but it is not very significant for the 
VRI bands. 

We have computed the reddening line for several of the most 
frequently used telescopes or instruments for KBO observations 


Table 1. Central wavelengths in A of different telescope filters 


Telescope 

B 

Eilters 

y R 

/ 

8.2m ESO-VLT 

4237 

5481 

6480 

7935 

3.6m ESO-NTT 

4212 

5442 

6416 

7950 

10m Keck 

4377 

5473 

6417 

7599 

8.2m Subaru 

4400 

5500 

6600 

8050 

3.6m CEHT 

4312 

5374 

6581 

8223 

Average 

4307.6 

5454 

6498.8 

8016.2 


and analyzed their differences. As seen in Fig.[2 the reddening 
line obtained for the Mould-BFR filters of the CEHT12k camera 
of the 3.6m CEHT telescope is higher, and the reddening lines 
obtained for the ESO EORSl camera of the 8.2m VET telescope 
or the ESO SuSi2 camera of the 3.6m NTT telescope are the 
lowest. The redenings for the LRIS camera on the 10m Keck 
telescope and the EOCAS cameras on the 8.2m Subaru telescope 
are intermediate. Therefore, interpretations based on a precise 
location of this line should be taken with caution. We have used 
the solar colors {B - y)o = 0.65, {V - R)q = 0.36, and (V - 
/)o = 0.70 from |Ramfrez et al.| ( |2012| l. Small changes in the 
solar colors do not affect the reddening line. 


5. Data analysis 


This work is focused on BVRI surface information or proper¬ 
ties of Centaurs and KBOs. We have analyzed the correlations 
between the color indexes B - V, V - R, R - I, V - I, B - I, 
and B - R, and the spectral gradient, Grt — also known as the 
slope parameter or reddening S , which expresses the re flectivity 
spectrum variation in percent of reddening per 1000 A (Hainaut 
|& Delsanti|2002] l — , and the correlations between these indexes 
and the absolute R-magnitude not corrected for phase effects, 
that is, = R(l, 1,0"). Correlations between these proper¬ 

ties and the orbital parameters perihelion, q, aphelion, Q, orbital 
inclination relative to the ecliptic, /, orbital eccentricity, e, and 
semi-major axis, a, have also been studied. For some families 
we have also analyzed the correlation with the Tisserand param¬ 
eter relative to Neptune, T v, the Tisserand parameter relative to 
Jupiter, Tj, the mean random impact speed velocity, Vrms, the or¬ 
bital excitation, e, and Opik’s if/. These parameters are explained 
below when used. As mentioned in Sect. we relied entirely on 
the published error bars, although for multiple measurements we 
checked for deviant points before carrying out the weighted av¬ 
eraging. Possible over- or underestimated errors for some data- 
points, most probable in the case of single measurements, are 
not expected to have any statistical influence on the methods we 
used and developed, unless if generalized. 

Each dynamical family, or group of families, was subject to 
84 tests for correlation. The exceptions were i) Centaurs, which 
were subjected to 92 tests due to the inclusion of Ty in the tests; 
ii) classical KBOs, which were subjected to 100 tests, given the 
inclusion of v„„s, A and and iii) resonant objects were only 
subjected to 76 tests, since correlations with T^i were not tested. 
These numbers are extremely important to correct for the false- 
discovery rate (EDR) as described in Sect. |3.3| 

We discuss the correlations for each dynamically family sep¬ 


arately from Sects. pTT| to 5.6 KBOs, as a whole, are discussed in 
Sect. mi binary or multiple objects in Sect. |5.8| and all objects 
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B-V 



Fig. 1. (S - V) vs. (V-R) — top — and (V - R) vs. (R-I) — bottom 
— reddening lines obtained with the most frequently used sets of BVRl 
filters in our database and the computed average reddening line. For 
the (B - V) vs. (V - R) case different filters have a strong effect on the 
location of the reddening line. For the (V-R) vs. (R -1) case the effect 
is almost negligible. 


— excluding Haumea family and retrogrades — are discussed in 
Sect. |5.9| For each case we also comment on the sample size and 
the consequent risk of missing correlations. Strictly speaking, 
this discussion is valid for a characterization of these families 
and groups under the assumption that the data are an unbiased 
sample of the parent population. The statistical tools, however, 
are general and can be applied to other classification systems. 


5.1. Centaurs 


to miss at a 10% risk p and at 2.5 cr is ppi.Sa- - 0.64. At the 
same 10% risk, but at 3 cr confidence level we have - 0.70. 
In other words, with the current sampling of Centaurs there is a 
risk higher than 10% of missing any possible strong correlation. 
Increasing the sampling of Centaurs is still very important. To 
have a risk lower than 10% of missing any correlation stronger 
than ai 0.5, we need a sample of n > 70 objects (more specifi¬ 
cally, pp 2 .Sa < 0.46 andp/j 3 o- < 0.50). 


5.2. Scattered disk objects 


The dynamical boundaries for the definition of scattered disk ob¬ 
jects, or SDOs, are ill defined. Different authors used slightly 
different definitions. Following the definition we used, we have 
n - 52 SDOs. Since one of them is a retrograde object and an¬ 
other one a Haumea family object, we have analyzed a total of 
n - 60 SDOs. Most color indexes correlate with each other (see 
bottom panels of Fig. [^and Table|^. The exception are (B - R) 
vs. (R - /) and (V - R) vs. (R - /), which indicates that it is 
difficult to extrapolate color values, although interpolating them 
seems possible. There is a correlation between Hr and q that we 
interpret as a biasing effect among SDOs, meaning that for orbits 
with high perihelia only the brightest objects can be detected, 
which was also discussed by Santos-Sanz et al.|(2012| l — for a 
review on biasing effects see |Kavelaars et al. ( 2()()8| l. At p = 0.1, 
with n — 61, we have pp 2 . 5 a- — 0.48 and pp^a- = 0.53, so some 
medium strength correlations might still be missed. A ~ 50% 
increase in the sampling of SDOs, to n = 90, would lower these 
values to pp 2 . 5 a- = 0.39 and pp^a- - 0.44. Note, however, that this 
biasing will not be removed solely by increasing the sample. 


5.3. Scattered or detached Kuiper Belt objects 

According to the classification scheme we have used, our sample 
has n = 7 detached KBOs (one of which is a binary object) and 
n = 15 scattered or detached KBOs (three of which are binaries). 
We opted to analyze them as one group (SDKBOs). The data 
show strong correlations between colors that we may consider 
as interpolations if we convert them into a spectrophotometric 
spectra. To illustrate, if we have a flat spectrum, the colors (B - 
R) or (V - /), for example, can ben seen as mere interpolations 
of the (B - I) color. Analogously, when in presence of a flat 
spectrum, the (B - /) color, for example, could be obtained by 
extrapolating from (B - R). This means that similarly as for the 
SDOs, consecutive colors or color extrapolations do not show 
evidence of correlation (see bottom panels of Fig. and Table 
1^. For example, (B - I) vs. (B - R) are strongly correlated, but 
the wavelength (B - R) is inside the (B - I) wavelength range, 
while (B - R) vs. (R - /) are not correlated. Because we are 
dealing with small-sampling values the p - 0.1 risk shows that 
it is easy to miss even very strong correlations among SDKBOs. 
For n = 20, pp2.5a- - 0.76 and pp^a- - 0.82. Clearly, it would be 
good to triple the sample. 


Our sample has n - 33 Centaurs. Because one of them was a 
retrograde object, however, we removed it and analyzed a total 
of n = 32 Centaurs. All of their surface colors indexes strongly 
correlate with each other, which indicates a strong degree of pre¬ 
dictability of one color when knowing another (see middle pan¬ 
els of Fig. 1^ and Table |^. No correlation between colors and 
absolute magnitude, and between these and any orbital parame¬ 
ter reaches a 2.5 cr level after correcting for FDR. We stress that 
with n - 32 objects, the lowest correlation value we expect not 
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5.4. Plutinos 

Our sample has n = 49 Plutinos, or (3 : 2) resonants, includ¬ 
ing Pluto and four binaries. All color indexes correlate with each 
other except for (V - R) vs. (R - I) (see Fig. and Table |^. 
Nonetheless, (B - V) correlated with (V - R), and (B - R) cor¬ 
relates with (V - I). Therefore, either color interpolation or ex¬ 
trapolations might be possible, but for the latter, Z?-band infor¬ 
mation would be required. No correlations between colors and 
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Fig. 2. (6 - V) vs. (V — K) plots — left column — and (V-R) vs. (R - /) plots — right column — including the reddening line. Solar colors 
are indicated with the O symbol. Top panels: All objects, Haumea family objects, and retrograde objects. Middle panels: Centaurs. Bottom panels: 
Scattered disk objects (SDOs), and scattered or detached Kuiper Belt objects (SDKBOs). The stronger correlation between (B - V) and (V — R) 
colors than between (V - R) and (R - /) is clear. Apart from the solar-like colors of Haumea family objects and the two color groups of Centaurs 
in the (B - V) vs. (V — R) there is no apparent relevant difference between the color distribution of each family of objects. 
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orbital parameters reach a 2.5 cr level after correcting for FDR. 
Given the non-detection of the (V-R)ws. (R-1) correlation with 
n - 42 objects, we highlight that at/3 = 0.1 risk p/? 2 . 5 cr = 0.56 
and Pysso- = 0.63. As for the SDOs, some medium-strength cor¬ 
relation might be missed here. Increasing the sample by ~ 50% 
is recommended. 


5.5. Other resonants 

Our sample has n - 12> non-Plutino resonants: n = 2 (1 : 1), 
or Neptune Trojans; n = 3 (5 : 4); n = 6 (4 : 3); n = 1 (11 : 
8 ); n = 1 (18 : 11); n = 10 (5 : 3) — 1 binary —; n = 1 
(12 : 7) — Haumea multiple system —; n - \ (19 : 11) — 
Haumea family member —; n - 15 (7 : 4) — 1 binary —; 
n = 1 (9 : 5); n = 10 (2 : 1) — 3 binaries —; n = 3 (9 : 4); 

n = 3 (7 : 3); n = 2 (12 : 5); « = 12 (5 : 2) — 1 binary —; 

n = 1 (8 : 3) — binary —; n - 2 {3 : 1); n = 1 (11 : 2); and 

n = 1 (11 : 3). We have removed the Haumea family members 

from the sample — (136108) Haumea and (019308) 1996T066 
— because of their extremely peculiar behaviors. This means 
that we analyzed the correlations for n - 11 objects (see Fig. 
[^and Table 1^. Absolutely all color indexes correlate with each 
other at a 3 cr confidence level after the FDR correction. The 
correlation magnitude weakens to p < 0.6 only for {B - R) vs. 
{R - I), (B - V) vs. (R - I), and (V - R) vs. {R - I). There is a 
moderate correlation between Hr and q, suggesting that we have 
an absolute magnitude bias among these objects similarly as for 
the SDOs. As for the previous cases, we detect no correlations 
between color indexes and orbital parameters. For n - 11, the 
sample size at which the lowest p is detected, the y 6 = 0.1 risk 
implies pp 2 . 5 a- = 0.44 and p/ 330 - = 0.50. Therefore, only weak 
or moderate correlations might be missed. As a whole, the other 
resonants do not to require a larger sampling to detect any highy 
relevant correlation. However, it is not likely that all the different 
resonances will possess the same surface color behaviors, and 
each resonance should be sampled such that they can be analyzed 
independently of the others. This still requires a large effort since 
even with n -20 objects one will only be able to be 90% certain 
to detect very strong correlations {pp 2 5 a- = 0.76 and p/ 330 - = 
0.81). 


5.6. Classical KBOs 


Classical KBOs, also known as classicals or CKBOs, are known 
for their color-inclination correlation, which is not clearly de¬ 
tected in any other family os KBOs (see Fig. top). As for the 
case of other resonants, we removed the seven Haumea fam¬ 
ily objects from the classicals sample: 2 OO 3 SQ 317 , 2003UZii7, 
2 OO 5 CB 79 , (024835) I 995 SM 55 , (055636) 2 OO 2 TX 300 , (120178) 
2 OO 3 OP 32 , and (145453) 2 OO 5 RR 43 . When plotting B - R vs. or¬ 
bital inclination i, we have three outliers with low B - R colors 
at very low inclination (see also Peixinho et al.|2008 1 . Roman- 
[ishin et ar] ( |2010| l analyzed a subfamily of classical objects with 
semi major-axes below the nominal location of the (3 : 2) res¬ 
onance, calling them inner classicals. That subfamily presented 
a bimodal behavior of surface colors, unlike the rest of classi¬ 
cal objects, and probably a distinct origin from the primordial 
protoplanetary disk. Although we did not use the same classifi¬ 
cation scheme, the outliers 1998WV24, and 2 OO 3 YL 179 belong 
to that group, and we decided to remove them from the analysis. 
There is, however, no similar indication why object 2002VDi3i 
might be an outlier., but because of its extreme outlier behavior 
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we chose to eliminate it as well. The final sample we analyzed 
has n = 117 classical objects (see Table]^. 


5.6.1. Correlation of color vs. color of classical KBOs 


CKBOs evidence neither (B-V) vs. (V-R) nor (V-R) vs. (R-I) 
correlations. And with their sample size, anything stronger than 
p ~ 0.43 should be detected at least at a 2.5 cr level, although 
within the p error-bars we see that these correlations among 
KBOs*, taken globally, increase when CKBOs are removed from 
the sample (see Table [^. The middle panel of Fig. shows that 
this result seems clear. Interestingly, there are more CKBOs in 
the {B - y) vs. (y - R) plot far below the reddening line than for 


any other family. We return to this question in Sect. 5.9.2 


5.6.2. Correlation of color vs. orbital parameters of classical 
KBOs 

The trend between the colors of classical objects and their orbital 
inclination was first reported by |Tegler & Romanishin ( 2000|l. 
Explicit computations of a correlation were made later by Tru- 
[jillo & Brown| ( [2002) l and |Hainaut & Delsanti| ( [2002| l. 

This correlation was thought to be caused by some colli- 
sional evolution mechanism, known as collisional resurfacing 
(CR) — see |Luu & Jewitt|fl996| ) — which seemed highly depen¬ 
dent on the orbital inclination because of its link with a stronger 
or weaker collisional environment. Different parameters have 
been used to analyze it, such as collisional velocity, orbital ex¬ 
citation, i/r parameter, or simply orbital inclination. We analyzed 
each one of them, also including the Tissserand parameter rela¬ 
tive to Neptune {Tff). 


Collisional velocity Vc'. The mean square of a KBO encounter 
velocity with its circular analog, that is, a non-inclined circular 
orbit with heliocentric distance a, is given by 




GMp 


(e2 + ^ 


V? 


(5) 


where e represents the orbital eccentricity, i stands for the or¬ 
bital inclination relative to the ecliptic, and a for the semi-major 
axis. If a is measured in AU, then vk ~ 29.8/ s/a [km For 
simplicity, we call this Vrms simply as collisional velocity Vc. 
Trujillo & Brown ( 2002|l discussed the possible link between 


the color-inclination correlation and impact velocity but, given 
the lack of detailed modeling available at the time, they also 
discussed the possibility that the correlation might be a con¬ 
sequence of multiple subpopulations with different primordial 


colors and inclinations. Stern (20021 explored the link between 
KBO colors and their mean collisional velocities in more de¬ 
tail and also accounted for the gravitational effect of each object 
using Vcoii - where Ug^c is the KBO’s escape ve¬ 

locity, calculated assuming a bulk density g - 1500 kgirP and a 
geometric albedo p = 0.04. Stem[ s Spearman-rank correlations 
of (B-R) vs. Vcoii and (B-V) vs. Vcoii werepBRv^^,, = -0.39 (p = 
0.005) and pbvucm - -0.44 (p - 0.0005), respectively, for a total 
of n = 81 KBOs. Note that these correlation values are weaker 


than those reported by Trujillo & Brown for colors vs. inclina¬ 
tion of classical KBOs also because Stern| analyzed all families 
together. Nonetheless, these works seemed to confirm the exis¬ 
tence of a physically meaningful correlation between KBO col¬ 
ors and their mean collision velocities, making a strong case for 
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Fig. 3. (6 - V) vs. (V - R) plots — left column — and {V — R) vs. (R - I) plots — right column — including the reddening line. Solar colors 
are indicated with the G symbol. Top panels: Plutinos and other resonants. Almost all non-Plutino resonants seem to lie on the reddening line in 
(B — V) vs. (V - R), but their stronger correlation value is within the error-bars of the correlations for other families. Middle panels: Classical 
KBOs and all other KBOs, except Haumea family and retrogrades. In contrast to all other families, CKBOs do not exhibit a(B - V) vs. (V — R) 
correlation. The existence of many CKBOs far below the reddening line (in BVR), which indicates the presence of many concave visible spectra, 
suggests that they may exhibit absorption features in the visible that only detailed spectroscopic studies could verify. Bottom panels: Binary or 
multiple KBOs. Article number, page 9 of 30 
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the CR scenario. From analyzing n = 50 CKBOs, Doressoundi- 


ram et al.|(|2002 |l repo rted Pbru„„ - -0.49 {p = 0.0004). Note 
s care in also accounting for Vesc to the 


that in spite of Stern 


collisional velocity from our sample, we measure a median dif¬ 
ference between Vcoii and Urms of only 0.02 km whereas the 
median Urms is 0.55 km This is a very small difference given 
the uncertainty on the albedos and densities, and is the reason 
why Vrms = Vc was used more often. It is noteworthy that Je- 
witt & Luu (2001|l, the proponents of the CR scenario, argued 


against it as the primary cause of color differences, since KBOs 
should then also exhibit large color variations with rotation, and 
observations indicated otherwise. iThebault & Doressoundiraml 
( [2003 | l and |Thebault|p003) found additional problems with this 
scenario. Through a numerical analysis, they computed the ki¬ 
netic energy KBOs would receive by collisions (hereafter Kcoii), 
considering several different swarms of impactors and also dif¬ 
ferent KBO distributions. Results showed correlations between 
Kcoii and perihelion, but also that Kcoii correlated systematically 
much better with eccentricity than with inclination. Observation- 
ally, there was evidence for color vs. perihelion and color vs. 
inclination correlations but, to date, no relevant correlations be¬ 
tween color and eccentricity has been found. Furthermore, from 
these simulations, Kcoii of Plutinos (3:2 resonants) was higher 
than the correlation of any other family, but observationally, their 
colors did not seem to differ. It is, nonetheless, true that this kind 
of numerical analysis has not been reanalyzed including the cur¬ 
rent understanding of the sculpting of the Kuiper Belt through 
migration processes. 


Orbital excitation e: [Hainaut & Delsanti] ( 2002| ) analyzed the 
color correlations with the object orbital excitation e, defined as 


e - Ve^ H- sin^ i. (6) 

This s parameter is related with Vrms (see Eq. [^, which is an¬ 
other estimate of the KBO velocity with respect to its circular 
analog. Using Pearson’s correlation coefficient (r^) instead of 
Spearman’s, they obtained for (B - R) vs. s of n=13 CKBOs 
r, = -0.77 (p = 0 .00:^, but no relevant values for any other 
families (r, < -0.1), nor for KBOs and Centaurs taken together 
(r, = -0.21, n=34). 


Opik’s if/: Santos-Sanz et al. ( 20091 analyzed the color correla- 
tions with a if/ parameter, taken from the theory of Opik |1976 i 
of close encounters of a small body with a planet, given by 


<A = 


5 0 • 2 • 

- + sin i 


(7) 


In our context, ^ relates with the average energy of the colli¬ 
sions experienced by a KBO. Again, the similarity between this 
parameter and the two previous ones is obvious. There is a sub¬ 
tle difference: -\JGMq \}/ is an approximation for the KBO aver¬ 
age encounter velocity with targets in circular orbits at the KBO 
instantaneous heliocentric distance on its elliptical orbit, that 
is, all local circular analogs. The two previous parameters con¬ 
sider only the impact with the global circular analog. Strangely, 
Santos-Sanz et al. ( 2009| l reported for {B - R) vs. if/ of n=64 CK¬ 


BOs p = -0.02. In our analysis we obtain P(b-r)iIi - -0-54);[j [Q 
(p < 10'®), as roughly expected given all the previous works. 


P-value is our estimate 
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Tisserand parameter T a?: The Tisserand parameter ( [Tisserand] 
1896| l, or invariant, has proven to be a useful tool in comet taxon¬ 
omy ( Kresak|1972[ Levison|1996 and references therein). It is a 
dynamical quantity kept approximately constant during a close 
encounter interaction between a planet and a minor body — for 
instance, a KBO could have its orbital parameters changed by a 
close interaction with Neptune, and yet the Tisserand parameter 
value would not change, allowing us to keep track of the same 
object. Since KBOs are mostly under gravitational influence of 
Neptune, we compute it relative to this planet using 


Tn - — + 2 (.Qg 

a y Uff 


( 8 ) 


where - 30.07 AU is the semi-axis of Neptune and a, e, 
and i are the orbital parameters of the KBO. A minor body on 
a circular orbit, coplanar with Neptune, with a - would 
have Tff - 3. A KBO, being a trans-Neptunian object, will have 
a > a^. If Tff > 3, then it will not cross the orbit of Neptune, 
whereas if < 3 means, generally, a Neptune-crosser. Usage of 
Tff in KBO taxonomy has been suggested by Elliot et al. ( 2005| l. 
Problems with the interpretation of Tff values arise when high in¬ 
clination values lead to a < 3, suggesting Neptune-crossing, 
even for impossible cases. Eor example, a KBO on a circular or¬ 
bit with a - 31 AU and i - 30° would have Tff - 2.7 and it 
would not cross the orbit of Neptune. Note that Tff parameter 
also relates with the relative encounter velocity between a KBO 
and Neptune: Vc oc yl3 — Tfj. 

Regardless of the dynamical implications, Jewitt et al. ( 2001) 
reported for a total of n = 15 CKBOs {B - R) vs. T ff and {U - B) 
vs. Tat correlation values of P(b-r)Th - -0.75;(Q2g (p - 0.0014), 
and P{u-b)Tn - -0-86j;Qj4 {p - 0.00004), respectively. 


Is there only one correlation? CKBOs show moderate to 
strong correlations, all of similar magnitudes, both with i, Vc, 
e, if/, and Tff, and a moderate to weak correlation with q. The 
confidence level for the correlation between color and eccentric¬ 
ity does not reach the 2.5 cr threshold after correcting for EDR, 
and no relevant correlation is detected between color and semi¬ 
major axis (see Tablej^. Erom the correlations = l.OOj;® ^®, 
p^^,e = 0.99;|;ggj, and p^y^ - 1.00;(qqj it is clear that both if/, s, 
and Vc will be correlated with any other variable with the same 
strength. Consequently, analyzing for correlations between col¬ 
ors and collisional velocity Vc, or orbital excitation e, or the Opik 
^ is equivalent and, from a statistical point of view, there is no 
reason to chose one over the other. Hereafter, we chose to work 
with Vc since it has been used more often in other works. 

Now, we consider the correlations between {B - R), Tff, Vc, 
i, and q more closely. Given the mutual dependence of most of 
these parameters, we have carried out a partial correlation anal¬ 
ysis to assess whether one of these correlation is dominant and 
the others are only spurious effects with no physical relevance, 
but are present because of their dependence on the dominant one 
(see Sect. [33- 

The partial correlations between {B - R), q, and i are 


I P(B-Rm = -0.41!“;[® (p - 0.00026; CL = 3.42 cr) 

1 P(B-R),.i = 0.14:°;;4 (p = 0.22; CL = 1.23 cr) 

showing that the correlation between {B - R) colors and orbital 
inclination is still present and significant after removing the ef¬ 
fect of the link between inclination and perihelion. Its strength. 
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Table 3. Spearman p correlation matrix between (6 - R) and orbital parameters for n = 19 CKBOs. 


B — R T]^ if/ 

Vc 

S 

i 

q 

e 

a 

B-R . 0.50!°;!^ -0-54!o;[o 

-0 54+0.12 


-0.11 

^•■^^-0.13 

-0.12 

-0 11+0-'4 

'+■ + +-0.14 

(p = 3 ■ 10-^) (p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p = 2 ■ 10-®) 

(p = 0.0011) 

(p = 0.0078) 

(p = 0.31) 

Tn . -0.78:°;°’ 

-0 77+0 07 
'+■' ' -0.05 


A oc+0.05 
^■°^-0.04 

0 

^•'^^-0.10 

^•^^-0.10 

^■^^-0.14 

(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p = 0.00083) 

(p = 0.026) 


1 

^•^^-0.01 

0 99+0.01 
'+■^^-0.01 

0 93+0.02 


0-68!°;°’ 

0 26+"°-'’ 
'+■^+’-0.13 


(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p = 0.034) 

Vc 

• 

1 

^•^^-0.01 

0 91+0.03 
^■^^-0.03 

-0.68!°;°° 

0 71+006 
' ^-0.07 

0.28!°;“ 



(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p = 0.0015) 

s 


• 

0 91+003 
'+'7' + -0.04 

-0.07 

A 79+0.06 

^ -0.07 

-0.12 




(p < 10-®) 

(p < 10-®) 

(p < 10-®) 

(p = 0.00052) 

i 



• 

_0 51+0.09 

A 47+0.08 

-0.09 

^■^■^- 0.12 





(p < 10-®) 

(p < 10-®) 

(p = 0.13) 

q 




• 

^•'^- 0.06 

A 04 + 0.12 
- 0.12 






(p < 10-®) 

(p = 0.69) 

e 





• 

^■^■^-0.09 
(p < 10-®) 

a 






• 


Notes. Correlation values represented as Significance p-values are two-tailed. 


however, diminishes to a more moderate value. Perihelion is not 
a function of inclination, and the fact that q and i are correlated 
could be due to some dynamical evolution mechanism in the way 
the Kuiper Belt was sculpted, or it may simply be the result of a 
sampling bias. 

The partial correlations between {B - R), T v, and q are 


I P(B-R)T,., = 0.39:°;!3 {p = 0.00035; CL = 3.57 tr) 

1 = 0.11!°;!4 {p = 0.35; CL = 0.93 tr) 

showing that the correlation between (B - R) colors and Tis- 
serand parameter is also still present and significant, although 
diminished to a more moderate value, after removing the effect 
of the link between Tn and perihelion, analogously to the previ¬ 
ous case. 

The partial correlations between {B - R), T v, and i are 


ip(B-RW- 0.17:0;[5 (p = 0.14; cl = 1.47 ct) 

1 P(B^R)i.T, = -0.17!°;!5 (p = 0.12; CL = 1.54fr) ’ 

which shows that none remains after eliminating the effect of the 
other. This should not be seen as strange given that by definition 
is highly dependent on /, presenting i=~o.s 5 , so it is very 
hard to separate one from the other. 

The partial correlations between {B - R), i, and Uc are 


I P(B^R)i.u, = (p = 0.59; CL = 0.53 a) 

1 P(B-R)v,.i = -0.20!°;|4 (p = 0.074; CL = 1.78 cr) ’ 

showing also that none remains after eliminating the effect of 
the other, even if the effect is stronger on inclination without the 


influence of collisional velocity than on Vc without the effect of 
i. Analogously to the previous case, by definition the collisional 
velocity depends on inclination {pvp=-o. 9 i), so it is hard to sepa¬ 
rate them. 

Finally, the partial correlations between {B - R), Tn, and Vc 
are 


I = -0.27!°;j4 (p = 0.016;CL = 2.41 cr) 

1 P(B-R)T,.v^ = (P = 0-14; CL = 1.48 cr) ’ 

indicating that none remains without the effect of the other vari¬ 
able, even if the case of p (^B-R)vc.Tfi is rather close to 2.5 cr. 

To conclude, for CKBOs 


- it is statistically equivalent to use i)/, e, or Vc- 

- the correlation between {B - R) and perihelia may be simply 
an effect that arises because in the current sample of mea¬ 
sured objects the perihelion is related with inclination, where 
inclination is the dominating variable. The same is true if we 
replace the variable inclination by the Tisserand parameter. 

- the correlations between {B - R) and both Tisserand parame¬ 
ter and inclination cannot be statistically separated from the 
close link between T^ and i. None dominates the other (see 

Fig- 01. 

- although less clear than for the previous case, the correlation 
between (B — R) and collisional velocity cannot be statisti¬ 
cally separated from the dependence of Vc on inclination and 
on Tn. 

- collisional velocity Vc, Tisserand parameter T^, or orbital 
inclination i may be the true physically relevant variables 
linked with (B - R) surface colors. 

- given previous studies on the lack of evidence that surface 
colors are affected by Vc, it is most likely that only either 
TN or i can be the variables that truly relate with colors — 


Article number, page 11 of 


30 







A&A proofs: manuscript no. peixinho+2015 



50 


40 


♦ Haumea family 
o KBOs* 

• Classical KBOs 



2.50 


♦ Haumea family 
o KBOs* 

• Classical KBOs 


0.80 1.20 1.60 2.00 2.40 

B-R 


Fig. 4. Top: (B - R) vs. orbital inclination plot of CKBOs and all KBOs 
with indication of Haumea family objects. The correlation for CKBOs 
alone is stronger than for all KBOs. The clustering of Haumea family 
objects had a strong effect on this correlation. Middle: Tisserand param¬ 
eter relative to Neptune (Tn) vs. orbital inclination as for the above. Tn 
is very strongly correlated with i for CKBOs. Bottom: (B - R] vs. Fiv 
as for the above. Given the very strong link between Tt^i and i, statis¬ 
tically one cannot know if only Tn or i have a (potentially) physically 
meaningful relation with (S - R) but none dominates the other. 


although causality cannot be assessed with this analysis —, 
and that dynamical history plays the central role. 


It is interesting to see that when the peculiar Haumea fam¬ 
ily objects are removed from the analysis of CKBOs, the color- 
inclination correlation diminishes considerably (see Fig. |^. 
|Peixinho et al.| ( |2008| l reported p - -0.70, for n - 69 objects, 
whereas now we have p - -0.50 for n - 19 objects, excluding 
the Haumea family. When Doressoundiram et al. ( 2005b| l ana¬ 
lyzed the correlation using two different classification schemes 
for CKBOs, one using q > 35 AU as cutoff and another using 
q > 39 AU, they reported p - -0.65 (n = 27) for the former 
and p - -0.41 (« = 20, and SL < 2cr) for the latter, four of the 
seven objects that had q <39 AU are now known to be members 
of the Haumea family. In other words, without the Haumea fam¬ 
ily, the color-inclination correlation among CKBOs is no longer 
a strong correlation but a moderate one. 

It is interesting to note that CKBOs do not show evidence 
for a correlation between {B - V) vs. (V - R), which is present 
in almost all the other more undersampled families. We return to 
this matter in Sect. 15.9.11 


5.7. All KBOs* 

We have also analyzed the n = 322 KBOs* in our sample, 
that is all objects except Centaurs, Haumea family objects, and 
retrograde orbits. All color indexes correlate with each other, 
as expected from previous works. However, it is interesting to 
note that the color-inclination correlation for all KBOs, although 
weaker, is within the error-bars of the same correlation for CK¬ 
BOs. This suggests that the colors of most non-CKBO families 
of KBOs are not totally uncorrelated with inclination, other¬ 
wise the correlation would be wiped out. The color-inclination 
correlation for n = 186 KBOs* removing the CKBOs is p = 
-0.30:!;jJ °’, p = 0.000042, CL = 4.09 cr (better than 3 cr after 
FDR correction). We must keep in mind that if the population’s 
p is -0.5, to have only a 10% risk of failing to detect it at a 2.5 cr 
level requires a sample of n = 55 objects, but if it is -0.3, we 
need n - 156. 

SDOs do have a sampling that should be enough to detect 
most correlations > | + 0.5|, but show P(B-R)i - “O-lO^^oJs 
(n - 54) with a p-value « 0.50. Plutinos have a smaller sam¬ 
pling and show a P{B-R)i — “0.30j;Q}5 (n = 45), where the first 
correlation is below the 2.5 cr level after correcting for FDR. For 
non-Plutino resonants we have P(B-R)i - “0-26tg[2 (n = 67), but 
below a 2 cr level after correcting for FDR. These results indicate 
that it is unlikely that any non-CKBO family will have a yet-to- 
be detected color-inclination correlation stronger than -0.5, but 
all may have a correlation of ~ -0.3 for which their sampling 
is too low to warrant detection. Clearly, the {B - R) vs. i trend 
is stronger among CKBOs, regardless of the current lower sam¬ 
pling of other families. 

To verify this behavior, we have simulated 1000 samples 
with 79 (ji:,y) data-points correlated with p = -0.50 (correspond¬ 
ing to our CKBOs) and 189 (.r,y) data-points randomly gener¬ 
ated, that is, with p ~ 0 (corresponding to the rest of our KBOs). 
The simulations show a distribution with {p)l‘^* - -0.17);° °^, 
which is not compatible with the observedp = -0.30j;gQ2, hence 
the weak color-inclination correlation seen among KBOs* is not 
simply a result of its stronger presence among CKBOs. 

Weak correlations of p ~ +0.2 are detected between col¬ 
ors and perihelion and eccentricity. Although weak in magnitude 
and only accounting for ~ 5% of the color variability — p^ may 
be roughly interpreted as the proportion of the variability of Y 
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that could be explained as a result of the variability of X — they 
are statistically significant, meaning that they are not a random 
effect. 

5.8. Binary or multiple KBOs 

Our sample has n - 48 binary or multiple objects, including 
Pluto and Haumea, and we have also analyzed the primaries of 
these systems as a group to see if they behave differently from 
the non-binaries. 

Their color-color correlation does not show a different be¬ 
havior from the previous cases. Their surface colors correlate 
both with Tn, with Uc, and i. A total of 26 out of these 48 objects 
are CKBOs, which are known to correlate with the aforemen¬ 
tioned variables. Only n - 35 binary or multiple KBOs possess 
measured (B - R) colors, which means « = 14 CKBOs. Analyz¬ 
ing only the n = 20 non-CKBOs, also excluding Haumea, we 
obtain p - -0.49;|;gjg, p - 0.027, CL - 2.21 cr (without correct¬ 
ing for FDR). The result is not significant, and given the sample 
size, we have ppi.Scr - -0.76, which means that the subsampling 
is too small to confirm the existence of P(B-R)i ~ -0.50, like for 
CKBOs, or P(B-R)i ~ -0.30 like for non-CKBO KBOs as dis¬ 
cussed above. 


5.9. Aif objects 

As indicated in Sect.j^ we have removed from our global anal¬ 
ysis objects belonging to the Haumea family and those with ret¬ 
rograde orbits. Our all* objects sample has n = 354 objects. The 
only relevant difference between the analysis of all* objects and 
that of KBOs* alone is the increase of the correlations between 
absolute magnitude Hr and both q and a. These two cases may 
be simply interpreted as a bias. That is, fainter objects, which are 
most often also smaller, are hard to observe at high semi-major 
axes or high perihelia. When we include Centaurs in the sam¬ 
ple, which are mostly small and faint, this becomes even more 
obvious. Nonetheless, it is a known fact that there are no large 
Centaurs, and this is not a bias. 


5.9.1. Consecutive color correlations 

Throughout all publications on the correlation between colors 
of Centaurs and KBOs, we see the tendency that with increas¬ 
ing wavelength of the region covered by the colors in question, 
the correlations tend to diminish. For example, in all works the 
correlations between (B - V) vs. {V - R) have been shown to 
be always stronger in magnitude than the correlations between 
i J - H) vs. (H - K). However, it is also a fact that the error-bars 
of visible colors are smaller than those of near-IR colors. Might 
this be the reason for the observed degradation of the correlation 
for increasing wavelength or not? 

To answer this question, we have analyzed the relative me¬ 
dian error-bar of each color normalized to the 6 cr range of the 
corresponding color. For example, the 6 cr range of{B-V) color 
for all* is 1.036, the median error-bar is 0.08, therefore, after nor¬ 
malization, if (B - V) ranges from 0 to 1 the median error-bars 
would be 0.077. With this procedure we can use Eq.|^and, with 
an assumption on the parent population p value, see the effect of 
the error-bars on the detectable p* (see Sect. |3.3[ ) 

Since the correlation of (B - V) vs. (V — R) usually is 
~ 0.5 - 0.6, we assume, for simplicity, two cases: i) the parent 
population p is 0.7 and it should be 0.7 regardless of the spec¬ 
tral region we examine, ii) the same as the previous, but with 


Table 4. Spearman p for consecutive colors. 




Observed 

Error 

Degraded 

Detectable 

Xvs. Y 

n 

pZ 

a-o 

Pq.i 

Fo.6 

Pl 32 . 5 cr 

All* 







(B-V)(V-R) 

292 

0.54!“;“ 

0.177 

0.54 

0.48 

0.22 

(V-R)(R-I) 

KBOs* 

280 

0 34+0-07 
-0.08 

0.192 

0.52 

0.47 

0.23 

{B-V)(V-R) 

266 

0 51+0-07 
^•'^^-0.07 

0.179 

0.53 

0.48 

0.23 

(V-R)(R-I) 

252 

-0.08 

0.195 

0.51 

0.47 

0.24 

KBOs* except CKBOs 





(B-V)(V-R) 

185 

0 61+0-07 
^•^^-0.08 

0.167 

0.55 

0.49 

0.28 

{V-R)(R-I) 

Classical KBOs 

175 

- 0.10 

0.179 

0.53 

0.48 

0.29 

(B-V)(V-R) 

78 

- 0.16 

0.234 

0.47 

0.43 

0.43 

{V-R)(R-I) 

76 

0.07!“;;“ 

0.241 

0.46 

0.42 

0.43 

Other resonants 






{B-V){V-R) 

67 

0.68!“;““ 

0.178 

0.54 

0.48 

0.50 

{V-R)(R-I) 

SDOs 

68 

0 44+0-13 

0.182 

0.53 

0.48 

0.56 

(B-V)(V-R) 

53 

^•^^- 0.14 

0.164 

0.55 

0.50 

0.51 

iV-R)(R-I) 

Plutinos 

45 

^•^^- 0.22 

0.175 

0.54 

0.49 

0.56 

(B-V)(V-R) 

45 

^•'^'^- 0.21 

0.186 

0.53 

0.47 

0.56 

{V-R)(R-I) 

42 

- 0.21 

0.265 

0.43 

0.40 

0.57 

Binary or multiple KBOs 





{B-V){V-R) 

34 

^•''^'^- 0.20 

0.147 

0.58 

0.51 

0.63 

{V-R)(R-I) 

Centaurs 

35 

^•*^-^- 0.23 

0.139 

0.59 

0.52 

0.65 

(B-V){V-R) 

26 

0 67+^-i^ 
^•^'-0.18 

0.119 

0.61 

0.54 

0.71 

(V-R)(R-I) 

SDKBOs 

28 

^•^^- 0.20 

0.159 

0.56 

0.50 

0.69 

{B-V){V-R) 

20 

^•'^'^- 0.40 

0.171 

0.55 

0.49 

0.70 

(V-R)(R-I) 

20 

^•^^- 0.35 

0.157 

0.56 

0.50 

0.70 


p - 0.6. Table [^summarizes the results. For each X vs. Y pair of 
colors we have their own sampling n, observed p, and normal¬ 
ized median error-bars cr^,. From Eq.j^we compute the degraded 
PDegmded = P* for both 0.7 and 0.6, and from the sampling n we 
solve Eq.j^to obtain the detectable p at a 2.5 cr level accepting 
only a riskyS = 0.1 of missing it, ppi.Scr (see Sect.|3.2|i. Statisti¬ 
cally p decreases with higher error-bars which are also typically 
associated with higher wavelengths but the variation is small. 
Observationally it decreases at a much faster rate. Therefore, we 
may conclude that, globally, the (B - V) vs. {V — R) correlation 
is stronger than the (V - R) vs. (R - 1) one, and it is not an effect 
created by larger error-bars, that is, the color-color correlations 
do diminish with increasing wavelength. This result does not dis¬ 
card the possibility that this cannot be the case for some specific 
subgroups of objects, however. 
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Examining these correlations for each family in more detail, 
we may see a similar behavior for all of them except for CKBOs. 
The sampling of Plutinos, binary or multiple KBOs, Centaurs, 
and SDKBOs does not warrant detection of correlations weaker 
than ~ 0.6 at a/3 risk of 10%. 

Since it is clear that the (V-R) vs. (R-I) correlation is indeed 
weaker than that of (B - V) vs. (V - R), it is not surprising that 
we do not detect it among SDOs, Plutinos, binary or multiple 
KBOs, or SDKBOs. For the latter we do not even detect the (B - 
V) vs. {V - R) correlation. For Centaurs, both are detected and 
show similar strengths. As a result of the low sampling, however, 
the correlation error-bars are too high to conclude that they are 
indeed similar, in contrast to what we have seen in other objects. 


The most interesting case, however, are CKBOs, since these 
objects show neither (B - V) vs. (V - R) nor (V - R) vs. (R - 
I) correlations (see Sect. 5.6. l| l. The visible spectrophotometric 
behavior of CKBOs, or at least one type of CKBOs, is distinct 
from other KBOs. We discuss it below. 


5.9.2. Peculiar convex spectrophotometric behaviors 


Surface spectra of Centaurs and KBOs tend to be straight in the 
visible and flatten toward the near-lR JHK bands, most partic¬ 
ularly when the visible slope is red, although the near-lR be¬ 
havior tends to be more erratic. Figures and show that in 
{B - V) vs. (V - R) the objects tend to lie in the vicinity of the 
reddening line, meaning that they possess rather straight spectra 
in the BVR wavelength range (pictorially; — or /), or above it, 
meaning that they possess a slightly convex spectra in that same 
range (pictorially: ^). In the (V - R) vs. {R - I) plots we may 
see that objects are distributed above and below the reddening 
line. We have identified 1/objects that fall far below the BVR 
reddening line: 2 OOOCF 104 , 2002 VTi 3 o, 2 OOIRZ 143 , 1999JDi32, 
2001HY65, 2 OOIHZ 58 , I 997 RT 5 , 2 OOOCF 105 , 2OO2VDi3i0(all 
Classical KBOs), 2 OOIKY 76 (Plutino), and 1999RJ215 (SDO) 
— these are all mostly classical KBOs. Their normalized spec¬ 
trophotometric reflectivity spectra are plotted in Fig. along 
with the spectra resulting from the mean colors of all* objects 
(the error bars from the latter are the standard deviation of the 
color distribution and not the average error on each color). Spec¬ 
tra are normalized to 1 at the center of the F-band Alter, and each 
one is shifted by 0.4 for clarity. The convex behavior these ob¬ 
jects present in BVR might be a mere consequence of a concave 
(pictorially: w) behavior in VRI, that is, the existence of an ab¬ 
sorption band at ~ 6000 - 7000 A. It is tempting to associate this 
behavior with aqueous altered minerals creating absorption fea¬ 
tures in this region, as once seen with spectroscopy on (38628) 
Huy a and (47932) 2 OOOGN 171 by |de Bergh erar]pn04| , but the 
effect on visible colors is not expected to be very strong (see 
|Barucci et al.||200^ for a review). The object 2 OOIRZ 143 does 
seem to be fully convex up to the I band, and we do not possess 
F-/ colors of objects 2002VDi3i, 2 OOOCF 104 , and 1999RJ215 to 
see if the spectral slope changes its convexity. Nonetheless, with 
the exception of 2 OOIRZ 143 , the other 10 aforementioned objects 
seem good candidates to measure for the long-sought hydration 
features among Centaurs and KBOs or any other absorption fea¬ 
ture in the visible spectrum. 


^ Note that 2002VD|3i was not included in the statistical analysis of 
classical KBOs because of its extreme outlier behavior. 
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Fig. 5. Spectrophotometric reflectivity of 1 1 peculiar objects exhibiting 
a convex behavior at ~ 4500 - 6000 A. Spectra are normalized to 1 
at 4 = 5454 A and are consecutively shifted by 0.4 for clarity. This 
apparent convexity may be created by a concave behavior at ~ 6000 - 
7000 A, which suggests some absorption feature on their surfaces in this 
wavelength range. The spectrum resulting from the mean B—V, V-R, 
and V-1 colors, and their standard deviations, for all* objects is shown 
in gray. 


5. 10. Bimodal colors of small objects 


IPeixinho et al.| ( |2012| l found that both Centaurs and KBOs 
present the so-called bimodal (B - R) visible color behavior at 
a p-value=0.001 (3.3 cr) for absolute magnitudes > 6.8. 

In a somewhat different manner, [Fraser & Brown] ( |2012] l found 
that all small objects (//*(«) > 6) with perihelion q < 35 AU 
are bimodal, which will therefore include Centaurs. With this 
new dataset we have carried out the same analysis as jPeixinho] 
et al. ( |2012[ ), obtaining a maximum of significance for the same 
bimodality for the n = 143 objects with Hufa) > 6.82, with 
a p-value=0.0046 (2.8 cr). Although the evidence for a bimodal 
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Fig. 6. (6 - R) colors vs. absolute magnitude (uncorrected for 

phase effects) for all objects in our database. Centaurs and KBOs are 
indicated with different symbols. Haumea family objects and the two 
retrograde objects are also highlighted. The Af-shape reported by |Peix-| 
|inho et aT] i 2012[ l is clearly visible, and faint object bimodality is de¬ 
tected for > 6.82, but at 2.8cr level (p-value=0.0046). 


behavior is very strong, we do not find a 3 cr result with the new 
dataset (see Fig.|^. 

Testing the n - 112 objects with perihelion q < 35 AU and 
H]i{a) > 6.0 from the all* sample results in a p-value=0.126 
(1.5 cr), which means that there is no evidence for bimodality. 

Interestingly, although the bimodality among Centaurs 
seems clear to the eye, when testing only for the n - 29 Cen¬ 
taurs in the sample we obtain a p-value=0.018 (2.4 cr). 


6. Conclusions 

We reanalyzed the correlations of BVRI-hand visible colors for 
n = 354 Centaurs and KBOs. The (B-V), (V-R), (R-I), (B-R), 
{V - I), and (B - I) color indexes were carefully computed from 
all the available measurements in the literature, discarding color 
indexes computed from two bands taken more than two hours 
apart, and rejecting clear outliers from the averaging process. 
The complete sample consisted of n = 366 objects, but we re¬ 
moved from the analysis n = 2 objects with retrograde orbits and 
n - 10 Haumea family objects. Nonetheless, these objects were 
also plotted in our figures. 

- We used an algorithm for the Spearman-rank correlation — 
also known as Spearman-p — taking into account not only 
the data error-bars but also the error on the estimate due to 
the finite sampling and corrected the significance levels using 
the false-discovery rate (FDR) due to the large number of 
tests each sample was being subjected to. 

- We analyzed for the first time the risk p of missing any pos¬ 
sible correlations as a result of our limited sampling and the 
sampling necessary to ensure detection at a chosen signifi¬ 
cance level. 

- With samples of n <70 objects, assumed to be unbiased 
extractions of the parent population, we have a /3 > 10% 
risk of missing strong correlations, and if n < 90 medium 
strength (moderate) correlations may still be missed. When 
analyzing the traditional dynamical families, only classical 


KBOs, as a whole, do not require larger sampling to detect 
presumably relevant correlations. If one aims to study the 
subfamilies of hot classical, cold classical, inner classicals, 
or other, then larger samples are still needed. Centaurs are 
particularly undersampled and should be more than doubled. 
Plutinos and SDOs would benefit from a 50% sample in¬ 
crease. Other resonants, as a whole, are not undersampled, 
but like the Plutinos, each resonant family should be ana¬ 
lyzed as a single group, and each one of them is largely un¬ 
dersampled. The sampling of the objects we grouped as scat¬ 
tered or detached KBOs would also need to be tripled. Note 
that the detection requirements for sample-size vs. correla¬ 
tion strength are general and would need to be applied to any 
possible subfamily or distinct classification one decides to 
study. However, the samples are assumed to be unbiased. In¬ 
cluding possible sample biasing effects would require a more 
complex analysis. 

- We studied the degradation of a detectable correlation as a 
function of the data error-bars and concluded that after a cer¬ 
tain correlation has been detected, it is pointless to increase 
the sampling while maintaining the same error-bars. Only by 
reducing the observational error-bars can a possibly stronger 
correlation be measured. 

- We showed that the location of so-called reddening line in 
the color-color plots is quite sensitive to the central wave¬ 
lengths of similar filters used on different telescopes. The ef¬ 
fect is strong in the (B - V) vs. {V - R) plots, but weak in the 
(V - R) vs. {R - /) plots. 

- The correlation between {B-V) and {V-R) colors is stronger 
than the correlation between {V - R) and {R - I), regardless 
of the observational error-bars. 

- We analyzed the different parameters that have been used 
in the context of the color-inclination correlation of classi¬ 
cal KBOs and the collisional resurfacing scenario using par¬ 
tial correlation tests. Statistically, it is equivalent to using 
the orbital excitation e, Opik’s i//, or the collisional velocity 
Vc- The detected color-perihelion correlation is most likely 
physically irrelevant and an effect of the strong relation be¬ 
tween perihelion q and inclination i among CKBOs. Both 
orbital inclination i and Tisserand parameter relative to Nep¬ 
tune Tat are the most likely physically relevant parameters 
related with the color variation of CKBOs, therefore the re¬ 
lation is probably related to dynamical history alone. Statis¬ 
tically, the collisional velocity cannot be separated from its 
dependence on inclination i or Tn- Without the Haumea fam¬ 
ily, the color-inclination correlation of CKBOs is no longer 
strong, but moderate {p ~ 0.5). 

- Although weaker than for CKBOs, a color-inclination corre¬ 
lation must still exist among the non-classical KBOs. 

- Classical KBOs are the only family that show neither {B-V) 
vs. {V - R) nor {V - R) vs. {R - I) correlation. 

- We identified ten objects, only two of which are not CKBOs, 
that exhibit a convex spectral behavior from their {B - V) and 
{V -R) colors and a concave behavior from their {V-R) and 
{R — I) colors, suggesting an absorption band or feature at 


~ 6000 - 7000 A, which renders them good candidates for 
detailed spectroscopic studies. 

We detected evidence for the {B-R) color bimodality for ob¬ 
jects with absolute magnitude HR{a) > 6.82 at a 2.8 cr level 
(p-value=0.0046) — see Peixinho et al. ( |2012| l. Paradoxi¬ 
cally, the evidence for bimodality among Centaurs reaches 
only 2.4 cr — see Peixinho et al. ( 2003| l —, and we find no 
evidence for bimodality of objects with perihelion q < 35 
AU and HR{a) ^ 6.0 — see [Fraser & Brown] ( |2012| l. 
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We must stress that we analyzed the correlations and sam¬ 
pling requirements of the traditional dynamical families of 
KBOs and Centaurs as defined in Sect.j^ with the aim to charac¬ 
terize each of them, and with the underlying statistical assump¬ 
tion that each sample is or will be unbiased. The physical and/or 
chemical relevance of the dynamical classification system used 
was not discussed. Using distinct classification systems, subsys¬ 
tems, or groupings should lead to distinct correlations as well 
as to distinct sampling considerations, and such works are to be 
encouraged. 

A similar study, extended to the near-IR colors, is under 
preparation. 
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— 3 cr FDR correction cutoff — 








(B-I) 

vs. 

(B-R) 

18 

0 79+0-12 

07378:?^ 

0.000085 

3.93 

2 

0.000326 

0.000071 

0.79 

0.84 

Grt 

vs. 

(V-D 

21 

0.000149 

3.79 

3 

0.000489 

0.000107 

0.75 

0.81 

Grt 

vs. 

(B-R) 

20 

^•'^-0.20 

-0.28 

0.000163 

3.77 

4 

0.000653 

0.000142 

0.76 

0.82 

Grt 

vs. 

(V-R) 

21 

0.000223 

3.69 

5 

0.000816 

0.000178 

0.75 

0.81 


continued on next page 
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Table 2. continued. 


FDR correction" 


x” 

VS. 

yc 

n"' 

(V-R) 

vs. 

(B-R) 

20 

(V-I) 

vs. 

(B-I) 

18 

(B-V) 

vs. 

(B-R) 

20 

(V-R) 

vs. 

(B-D 

18 

— 2.5 cr FDR correction cutoff 

Plutinos 




Grt 

vs. 

(B-D 

18 

(B-I) 

vs. 

(B-R) 

18 

Grt 

vs. 

(V-I) 

21 

Grt 

vs. 

(B-R) 

20 

Grt 

vs. 

(V-R) 

21 

(V-R) 

vs. 

(B-R) 

20 

(V-I) 

vs. 

(B-I) 

18 

(B-V) 

vs. 

(B-R) 

20 

(V-R) 

vs. 

(B-I) 

18 

(B-V) 

vs. 

(B-I) 

18 

(V-I) 

vs. 

(B-R) 

20 

(V-R) 

vs. 

(V-I) 

21 

(R-I) 

vs. 

(V-I) 

21 

Grt 

vs. 

(B-V) 

20 

(B-V) 

vs. 

(V-R) 

20 

(B-R) 

vs. 

Vc 

20 

(B-V) 

vs. 

(V-D 

20 

— 3 cr FDR correction cutoff — 

(V-I) 

vs. 

Vc 

22 

(R-I) 

vs. 

(B-I) 

18 

(V-I) 

vs. 

i 

22 

— 2.5 cr FDR correction cutoff 

Other resonants 


B-I 

vs. 

B-R 

62 

(B-V) 

vs. 

(B-I) 

62 

(B-V) 

vs. 

(B-R) 

67 

(B-V) 

vs. 

(V-I) 

64 

(B-V) 

vs. 

(V-R) 

67 

Grt 

vs. 

(B-I) 

62 

Grt 

vs. 

(B-R) 

67 

Grt 

vs. 

(B-V) 

67 

Grt 

vs. 

(R-I) 

68 

Grt 

vs. 

(V-D 

68 

Grt 

vs. 

(V-R) 

71 

(R-I) 

vs. 

(B-D 

62 

(R-I) 

vs. 

(V-D 

68 

(V-I) 

vs. 

(B-I) 

62 

(V-I) 

vs. 

(B-R) 

64 

(V-R) 

vs. 

(B-D 

62 

(V-R) 

vs. 

(B-R) 

67 

(V-R) 

vs. 

(V-D 

68 

(R-I) 

vs. 

(B-R) 

64 

(B-V) 

vs. 

(R-I) 

64 

(V-R) 

vs. 

(R-I) 

68 




0 73 + 0.16 
'^■'■’-0.32 

'-^-0.26 
0 71 + 0-'7 

07 781 

^■^^-0.36 


0 . 82!|;;^2 

0 79 + 0-12 

0 7 378:?^ 

07578:?^ 

^■'^-0.20 

-0.28 
0 73+o TS 
^■^*^-0.32 

-0.26 
0 71+^1^ 

07 781 
067^8:^^ 

'7-^'-Q.3Q 

0 62+°'* 
-Q.27 

0.58+ 

0 57+° 

781 


78:1 


0.56^ 

0.55^ 

-0.54^ 

0.5F 


:8:i§ 

0.28 


-0 45 + 0-28 
^•'^^-0.33 

-0 45+^-^^ 
^■^'^-0.21 


0.86!8;°° 

0.80 
0.69_, 
0 . 68 +' 
0 . 86 _ 
0.83+ 
0.75+ 
0.61+ 


8:At 

8:A8 

8:4? 


ai 

J.14 
0 79+0.07 

07^+8:4§ 

070^8:44 

-Q-iQ 

0.65!;}° 

0-83!:| 

0 75+°-0» 

: 78 : 4 § 


0.75+ 

0.76;^ 

0 . 66 + 


I.IO 

78:4a 

0 55+8''8 

0 53+°'° 

, 78 : 1 ? 


0.44' 


— 3 cr FDR correction cutoff — 

Hr vs. q 71 

— 2.5 cr FDR correction cutoff — 


0.15 


-0 42+°'° 

-0.09 


Classical KBOs 

continued on next page 


P-value^ CL{o-y 

0.000257 3.65 

0.000352 3.57 

0.000396 3.54 

0.001018 3.29 


0.000033 

4.15 

0.000085 

3.93 

0.000149 

3.79 

0.000163 

3.77 

0.000223 

3.69 

0.000257 

3.65 

0.000352 

3.57 

0.000396 

3.54 

0.001018 

3.29 

0.002264 

3.05 

0.003723 

2.90 

0.006273 

2.73 

0.006444 

2.72 

0.009851 

2.58 

0.011730 

2.52 

0.014885 

2.44 

0.023125 

2.27 

0.035165 

2.11 

0.035576 

2.10 

0.035680 

2.10 


< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

< 10-° 

> 8.0 

0.000003 

4.70 

0.000006 

4.53 

0.000197 

3.72 

0.000246 

3.67 


D 

i 2.5 cr' 

m 

I-3crj 

m 

P02.5a-‘ 

6 

0.000979 

0.000213 

0.76 

7 

0.001142 

0.000249 

0.79 

8 

0.001305 

0.000284 

0.76 

9 

0.001468 

0.000320 

0.79 


1 

0.000163 

2 

0.000326 

3 

0.000489 

4 

0.000653 

5 

0.000816 

6 

0.000979 

7 

0.001142 

8 

0.001305 

9 

0.001468 

10 

0.001632 

11 

0.001795 

12 

0.001958 

13 

0.002121 

14 

0.002284 

15 

0.002447 

16 

0.002611 

17 

0.002774 


0.000036 

0.79 

0.000071 

0.79 

0.000107 

0.75 

0.000142 

0.76 

0.000178 

0.75 

0.000213 

0.76 

0.000249 

0.79 

0.000284 

0.76 

0.000320 

0.79 

0.000355 

0.79 

0.000391 

0.76 

0.000426 

0.75 

0.000462 

0.75 

0.000497 

0.76 

0.000533 

0.76 

0.000568 

0.76 

0.000604 

0.76 


18 0.002937 

19 0.003100 

20 0.003263 


0.000639 0.74 

0.000675 0.79 

0.000711 0.74 


1 

0.000163 

0.000036 

0.47 

2 

0.000326 

0.000071 

0.47 

3 

0.000489 

0.000107 

0.45 

4 

0.000653 

0.000142 

0.46 

5 

0.000816 

0.000178 

0.45 

6 

0.000979 

0.000213 

0.47 

7 

0.001142 

0.000249 

0.45 

8 

0.001305 

0.000284 

0.45 

9 

0.001468 

0.000320 

0.45 

10 

0.001632 

0.000355 

0.45 

11 

0.001795 

0.000391 

0.44 

12 

0.001958 

0.000426 

0.47 

13 

0.002121 

0.000462 

0.45 

14 

0.002284 

0.000497 

0.47 

15 

0.002447 

0.000533 

0.46 

16 

0.002611 

0.000568 

0.47 

17 

0.002774 

0.000604 

0.45 

18 

0.002937 

0.000639 

0.45 

19 

0.003100 

0.000675 

0.46 

20 

0.003263 

0.000711 

0.46 

21 

0.003426 

0.000746 

0.45 

22 

0.003589 

0.000782 

0.44 


PyS3o-* 

0.82 

0.84 

0.82 

0.84 


0.84 

0.84 

0.81 

0.82 

0.81 

0.82 

0.84 

0.82 

0.84 

0.84 

0.82 

0.81 

0.81 

0.82 

0.82 

0.82 

0.82 

0.80 

0.84 

0.80 


0.53 

0.53 

0.51 

0.52 

0.51 

0.53 

0.51 

0.51 

0.51 

0.51 

0.50 

0.53 

0.51 

0.53 

0.52 

0.53 

0.51 

0.51 

0.52 

0.52 

0.51 

0.50 
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Table 2. continued. 


FDR correction" 


x” 

VS. 

yc 

n"' 

(B-R) 

vs. 

s 

79 

(B-R) 

vs. 


79 

(B-R) 

vs. 

Vc 

79 

(B-V) 

vs. 

(B-R) 

78 

Grt 

vs. 

(V-R) 

94 

(B-I) 

vs. 

(B-R) 

53 

(B-R) 

vs. 

i 

79 

(B-R) 

vs. 

Tn 

79 

Grt 

vs. 

(V-I) 

80 

Grt 

vs. 

(B-R) 

78 

(B-V) 

vs. 

s 

78 

(B-V) 

vs. 

Vc 

78 

(B-V) 

vs. 


78 

(V-R) 

vs. 

Tn 

97 

(B-V) 

vs. 

(B-I) 

53 

Grt 

vs. 

Vc 

94 

Grt 

vs. 

Tn 

94 

Grt 

vs. 

-A 

94 

(V-R) 

vs. 

(B-R) 

78 

Grt 

vs. 

s 

94 

Grt 

vs. 

(B-I) 

53 

(V-R) 

vs. 

-A 

97 

(B-V) 

vs. 

i 

78 

(V-R) 

vs. 

Vc 

97 

(V-R) 

vs. 

i 

97 

(V-R) 

vs. 

s 

97 

(B-I) 

vs. 

Vc 

53 

(B-I) 

vs. 

s 

53 

Grt 

vs. 

i 

94 

(B-I) 

vs. 

-A 

53 

Hr 

vs. 

Tn 

98 

(V-I) 

vs. 

(B-I) 

53 

— 3 cr FDR correction cutoff — 

(V-R) 

vs. 

(V-I) 

81 

(B-R) 

vs. 

q 

79 

(B-V) 

vs. 

Tn 

78 

Hr 

vs. 

i 

98 

Hr 

vs. 

-A 

98 

Hr 

vs. 

Vc 

98 

Hr 

vs. 

s 

98 

(B-I) 

vs. 

i 

53 

(V-I) 

vs. 

(B-R) 

65 

(R-I) 

vs. 

(V-D 

76 

— 2.5cr 

FDR correction cutoff 

Binary or multiple KBOs 


(B-I) 

vs. 

(B-R) 

28 

(B-V) 

vs. 

(B-R) 

34 

Grt 

vs. 

(B-R) 

34 

Grt 

vs. 

(V-D 

37 

Grt 

vs. 

(V-R) 

38 

(V-R) 

vs. 

(B-R) 

34 

(V-I) 

vs. 

(B-R) 

34 

(B-V) 

vs. 

(B-I) 

28 

Grt 

vs. 

(B-I) 

28 


P-7 

P-value^ 

CL(cry 

D 

i 2.5 cr' 

m 

I-3cri 

m 

Ppi.Sc/ 

PyS3o-' 

-0.53:“-|+ 

< 10-“ 

5.06 

1 

0.000124 

0.000027 

0.42 

0.47 

-0 54+“-'2 

< 10-“ 

5.12 

2 

0.000248 

0.000054 

0.42 

0.47 

-0 54+^-^^ 

10 

< 10-“ 

5.12 

3 

0.000372 

0.000081 

0.42 

0.47 

15 

< 10-“ 

5.22 

4 

0.000496 

0.000108 

0.42 

0.47 

0 54+012 

< 10“ 

> 8.0 

5 

0.000620 

0.000135 

0.39 

0.43 

0.60!-; 

0.000002 

4.78 

6 

0.000744 

0.000162 

0.51 

0.57 

-0.50!-! 

0.000002 

4.73 

7 

0.000868 

0.000189 

0.42 

0.47 

0.50+I 

0.000003 

4.71 

8 

0.000992 

0.000216 

0.42 

0.47 

0.48:- 

0.000006 

4.53 

9 

0.001116 

0.000243 

0.42 

0.47 

0.47+- 2 

0.000017 

4.30 

10 

0.001240 

0.000270 

0.42 

0.47 

-0 45+^-'^ 

0.000041 

4.10 

11 

0.001364 

0.000297 

0.42 

0.47 

-0.44+“-[2 

0.000052 

4.05 

12 

0.001488 

0.000324 

0.42 

0.47 

-0.44+i ^ 

0.000063 

4.00 

13 

0.001612 

0.000351 

0.42 

0.47 

0.39! + 

0.000084 

3.93 

14 

0.001736 

0.000378 

0.38 

0.43 

0-51!! 

0.000104 

3.88 

15 

0.001860 

0.000405 

0.51 

0.57 

-0.38!-[ 

0.000159 

3.78 

16 

0.001984 

0.000432 

0.39 

0.43 

0.38+|j“ 

0.000167 

3.76 

17 

0.002108 

0.000459 

0.39 

0.43 

-0.38! + 

0.000168 

3.76 

18 

0.002232 

0.000486 

0.39 

0.43 

0.41+“-[^ 

0.000186 

3.74 

19 

0.002356 

0.000513 

0.42 

0.47 

-0.38! + 

0.000187 

3.74 

20 

0.002480 

0.000540 

0.39 

0.43 

0.49!“- 

0.000205 

3.71 

21 

0.002604 

0.000567 

0.51 

0.57 

-0 36+■ 

0.000264 

3.65 

22 

0.002728 

0.000594 

0.38 

0.43 

-0.40!-| 

0.000265 

3.65 

23 

0.002852 

0.000621 

0.42 

0.47 

-0.36+“- 2 

0.000303 

3.61 

24 

0.002976 

0.000648 

0.38 

0.43 

-0.36^ 

0.000323 

3.60 

25 

0.003100 

0.000675 

0.38 

0.43 

-0.36^1! 

0.000356 

3.57 

26 

0.003224 

0.000702 

0.38 

0.43 

-0.47!-| 

0.000357 

3.57 

27 

0.003348 

0.000729 

0.51 

0.57 

-0.47!“-1 

0.000359 

3.57 

28 

0.003472 

0.000756 

0.51 

0.57 

-0 36+ ' 

0.000390 

3.55 

29 

0.003596 

0.000783 

0.39 

0.43 

-0.47!“- 

0.000417 

3.53 

30 

0.003720 

0.000810 

0.51 

0.57 

0 34+0-09 

0.000507 

3.48 

31 

0.003844 

0.000837 

0.38 

0.43 

0.45i| 

0.000834 

3.34 

32 

0.003968 

0.000864 

0.51 

0.57 

0.36!“-!3 

0.000940 

3.31 

33 

0.004092 

0.000891 

0.42 

0.47 

0 36+“ '+ 

0.001117 

3.26 

34 

0.004216 

0.000918 

0.42 

0.47 

0 36+“- 

u.JU_o 14 

0.001192 

3.24 

35 

0.004340 

0.000945 

0.42 

0.47 

-0.32!“';“ 

0.001426 

3.19 

36 

0.004464 

0.000972 

0.38 

0.43 

-0.32+i'^ 

0.001585 

3.16 

37 

0.004588 

0.000999 

0.38 

0.43 


0.001940 

3.10 

38 

0.004712 

0.001026 

0.38 

0.43 

-0 30+ ■' 

0.002265 

3.05 

39 

0.004836 

0.001053 

0.38 

0.43 

-0.41!“'[° 

0.002328 

3.04 

40 

0.004960 

0.001080 

0.51 

0.57 

0-36!i[^ 

0.003186 

2.95 

41 

0.005084 

0.001107 

0.46 

0.52 

0 33+^+^ 
u.JJ_o.l7 

0.004019 

2.88 

42 

0.005208 

0.001134 

0.43 

0.48 


continued on next page 


0.78!;! 

0.74!;[f 

0 76 ^“'° 
-Qifi 

0 76 +“ '° 
-Q15 

0 72 +“ '^ 
-0.18 

0 70+“ '' 
'+■''+-0.16 

0 75+“ '' 
'+■'-+-0.18 

0 75+“ '' 

'+■'-+-0.17 


< 10-“ 
< 10-“ 
< 10 -“ 
< 10 -“ 
< 10 -“ 
0.000002 
0.000004 
0.000005 
0.000005 


5.42 
> 8.0 
5.06 
> 8.0 
> 8.0 
4.78 
4.60 
4.58 
4.57 


1 

2 

3 

4 

5 

6 
7 


0.000148 

0.000295 

0.000443 

0.000590 

0.000738 

0.000886 

0.001033 

0.001181 

0.001329 


0.000032 

0.000064 

0.000096 

0.000129 

0.000161 

0.000193 

0.000225 

0.000257 

0.000289 


0.67 

0.62 

0.62 

0.60 

0.59 

0.62 

0.62 

0.67 

0.67 


0.73 

0.68 

0.68 

0.66 

0.65 

0.68 

0.68 

0.73 

0.73 
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Table 2. continued. 


FDR correction" 


x” 

VS. 

yc 

n"' 

(V-R) 

vs. 

(V-I) 

37 

(V-I) 

vs. 

(B-D 

28 

(V-R) 

vs. 

(B-I) 

28 

(B-R) 

vs. 

In 

35 

(B-R) 

vs. 

Vc 

35 

Grt 

vs. 

Tn 

38 

(B-R) 

vs. 

i 

35 

(B-I) 

vs. 

Vc 

28 

(B-V) 

vs. 

(V-D 

33 

Grt 

vs. 

(B-V) 

34 

(V-I) 

vs. 

Tn 

47 

— 3 cr FDR correction cutoff — 

(V-R) 

vs. 

Tn 

38 

(B-V) 

vs. 

Vc 

34 

(B-V) 

vs. 

(V-R) 

34 

(B-I) 

vs. 

Tn 

28 

Grt 

vs. 

Vc 

38 

(B-V) 

vs. 

Tn 

34 

Grt 

vs. 

(R-D 

35 

Hr 

vs. 

i 

38 

(V-R) 

vs. 

Vc 

38 

(V-I) 

vs. 

Vc 

47 

(R-D 

vs. 

(V-I) 

35 

Grt 

vs. 

i 

38 

(V-R) 

vs. 

i 

38 

(B-V) 

vs. 

i 

34 

(B-I) 

vs. 

i 

28 

— 2.5 cr FDR correction cutoff 

KBO* 




(B-I) 

vs. 

(B-R) 

207 

(B-I) 

vs. 

Tn 

209 

(B-I) 

vs. 

Vc 

209 

(B-R) 

vs. 

Tn 

268 

(B-R) 

vs. 

Vc 

268 

(B-R) 

vs. 

i 

268 

(B-V) 

vs. 

(B-I) 

209 

(B-V) 

vs. 

(B-R) 

264 

(B-V) 

vs. 

(R-D 

223 

(B-V) 

vs. 

Tn 

266 

(B-V) 

vs. 

(V-D 

233 

(B-V) 

vs. 

(V-R) 

266 

(B-V) 

vs. 

Vc 

266 

(B-V) 

vs. 

i 

266 

Grt 

vs. 

(B-D 

209 

Grt 

vs. 

(B-R) 

264 

Grt 

vs. 

(B-V) 

266 

Grt 

vs. 

(R-D 

251 

Grt 

vs. 

Tn 

294 

Grt 

vs. 

(V-D 

261 

Grt 

vs. 

(V-R) 

294 

Grt 

vs. 

Vc 

294 

Grt 

vs. 

i 

294 

Hr 

vs. 

q 

291 


P-7 

P-value^ 

CL(cry 

D 

i 2.5 cr' 

m 

i3 + 

m 

PI32.s7 

PyS3o-' 

o.67+;;;^ 

0.000005 

4.56 

10 

0.001476 

0.000321 

0.60 

0.66 

0.69!|'^ 

0.000045 

4.08 

11 

0.001624 

0.000354 

0.67 

0.73 

0.69!i| 

0.000050 

4.05 

12 

0.001771 

0.000386 

0.67 

0.73 

0 62+^-^ 

0.000063 

4.00 

13 

0.001919 

0.000418 

0.61 

0.68 

-0.62+^- 1 

0.000067 

3.98 

14 

0.002067 

0.000450 

0.61 

0.68 

15 

0.000245 

3.67 

15 

0.002214 

0.000482 

0.59 

0.65 

-0.58![5 

0.000260 

3.65 

16 

0.002362 

0.000514 

0.61 

0.68 

-0.63!“[° 

0.000306 

3.61 

17 

0.002510 

0.000546 

0.67 

0.73 

0 59+“-^ 

0.000331 

3.59 

18 

0.002657 

0.000579 

0.63 

0.69 

0 57+^-^^ 

C'--’ '_n JQ 

0.000410 

3.53 

19 

0.002805 

0.000611 

0.62 

0.68 

0 49+0-13 

0.000446 

3.51 

20 

0.002952 

0.000643 

0.54 

0.60 

0.53!°[3 

0.000706 

3.39 

21 

0.003100 

0.000675 

0.59 

0.65 

-0 55^ ■' 

0.000750 

3.37 

22 

0.003248 

0.000707 

0.62 

0.68 

0.55i'5 

0.000815 

3.35 

23 

0.003395 

0.000739 

0.62 

0.68 

0 5978:?§ 

0.000858 

3.33 

24 

0.003543 

0.000771 

0.67 

0.73 

_0 51+0.15 

0.001201 

3.24 

25 

0.003690 

0.000804 

0.59 

0.65 

0 52+^-^^ 

14 

0.001522 

3.17 

26 

0.003838 

0.000836 

0.62 

0.68 

0.50!!’ 

0.002226 

3.06 

27 

0.003986 

0.000868 

0.61 

0.68 

-0.48+11 

0.002353 

3.04 

28 

0.004133 

0.000900 

0.59 

0.65 

-0.48+°- ' 

0.002437 

3.03 

29 

0.004281 

0.000932 

0.59 

0.65 

-0.43!|| 

0.002467 

3.03 

30 

0.004429 

0.000964 

0.54 

0.60 

A 40 + 0:20 
27 

0.002609 

3.01 

31 

0.004576 

0.000996 

0.61 

0.68 

-0.47!°[® 

0.003044 

2.96 

32 

0.004724 

0.001029 

0.59 

0.65 

-0 46+ ■' 

0.003376 

2.93 

33 

0.004871 

0.001061 

0.59 

0.65 

-0.48!°;! 

0.003660 

2.91 

34 

0.005019 

0.001093 

0.62 

0.68 

-0 53+^-i^ 

^•^‘^-0.15 

0.003746 

2.90 

35 

0.005167 

0.001125 

0.67 

0.73 


0.83!°' 

0.36! 

-0.39! 

0.41! 

-0.42! 

-0.39! 

0.75! 

0.76! 

0.41! 

0.35! 

0.58! 

0.51! 

-0.37! 

-0.34! 

0.77! 

0.74! 

0.59! 

0.48! 

0.36! 

0.72! 

0.71! 

-0.36! 

-0.35! 

-0.37!;, 


< 10-® 
< 10-® 
< 10-® 
< 10 -® 
< 10-® 
< 10 -® 
< 10 -® 
< 10 -® 
< 10-® 
< 10-® 
< 10-® 
< 10-® 
< 10 -® 
< 10 ^® 
< 10 -® 
< 10 -® 
< 10-® 
< 10^® 
< 10-® 
< 10-® 
< 10-® 
< 10-® 
< 10-® 
< 10 -® 


> 8.0 
5.42 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 
> 8.0 


1 

2 

3 

4 

5 

6 

7 

8 
9 


0 
0 
0 , 
0 
0 
0 
0 
0 
0 , 

10 0 , 
11 0 , 
12 0 , 

13 0, 

14 0, 

15 0 

16 0 

17 0, 

18 0, 

19 0, 

20 0 .' 
21 0 .' 
22 0 

23 0 

24 0.' 


000148 

000295 

000443 

000590 

000738 

000886 

001033 

001181 

001329 

001476 

001624 

001771 

001919 

002067 

002214 

002362 

002510 

002657 

002805 

002952 

003100 

003248 

003395 

003543 


0.000032 

0.000064 

0.000096 

0.000129 

0.000161 

0.000193 

0.000225 

0.000257 

0.000289 

0.000321 

0.000354 

0.000386 

0.000418 

0.000450 

0.000482 

0.000514 

0.000546 

0.000579 

0.000611 

0.000643 

0.000675 

0.000707 

0.000739 

0.000771 


0.26 

0.26 

0.26 

0.23 

0.23 

0.23 

0.26 

0.23 

0.25 

0.23 

0.25 

0.23 

0.23 

0.23 

0.26 

0.23 

0.23 

0.24 

0.22 

0.23 

0.22 

0.22 

0.22 

0.22 


0.30 

0.29 

0.29 

0.26 

0.26 

0.26 

0.29 

0.26 

0.29 

0.26 

0.28 

0.26 

0.26 

0.26 

0.29 

0.26 

0.26 

0.27 

0.25 

0.26 

0.25 

0.25 

0.25 

0.25 
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Table 2. continued. 


x” 

VS. 

yc 

n"' 

F 

(R-I) 

vs. 

(B-D 

206 

0.56 

(R-I) 

vs. 

(B-R) 

222 

0.44 

(R-I) 

vs. 

(V-D 

253 

0.53 

(V-I) 

vs. 

(B-D 

209 

0.73 

(V-D 

vs. 

(B-R) 

232 

0.68 

(V-I) 

vs. 

Tn 

281 

0.33 

(V-D 

vs. 

Vc 

281 

-0.34 

(V-D 

vs. 

i 

281 

-0.31 

(V-R) 

vs. 

(B-D 

209 

0.65 

(V-R) 

vs. 

(B-R) 

264 

0.67 

(V-R) 

vs. 

Tn 

297 

0.32 

(V-R) 

vs. 

(V-D 

262 

0.59 

(V-R) 

vs. 

Vc 

297 

-0.32 

(V-R) 

vs. 

i 

297 

-0.32 

(B-D 

vs. 

i 

209 

-0.33 

(V-R) 

vs. 

(R-D 

252 

0.30 

(B-R) 

vs. 

e 

268 

-0.25 

(B-R) 

vs. 

q 

268 

0.25 

Grt 

vs. 

q 

294 

0.22 

Grt 

vs. 

e 

294 

-0.21 

(B-V) 

vs. 

q 

266 

0.22 

(B-I) 

vs. 

q 

209 

0.24 

(B-V) 

vs. 

e 

266 

-0.21 

(B-I) 

vs. 

e 

209 

-0.24 

(V-D 

vs. 

e 

281 

-0.21 

(V-I) 

vs. 

q 

281 

0.20 

— 3 cr FDR correction cutoff — 


(R-I) 

vs. 

Vc 

253 

-0.20 

(R-I) 

vs. 

Tn 

253 

0.19 

Hr 

vs. 

a 

291 

-0.18 

(V-R) 

vs. 

q 

297 

0.17 

(V-R) 

vs. 

e 

297 

-0.16 

Hr 

vs. 

e 

291 

0.16 

(R-I) 

vs. 

q 

253 

0.17 

— 2.5 cr FDR correction cutoff - 

— 

Air 





(B-D 

vs. 

(B-R) 

226 

0.84 

(B-D 

vs. 

Vc 

228 

-0.36 

(B-D 

vs. 

i 

228 

-0.33 

(B-R) 

vs. 

Tn 

297 

0.29 

(B-R) 

vs. 

Vc 

297 

-0.41 

(B-R) 

vs. 

i 

297 

-0.39 

(B-V) 

vs. 

(B-D 

228 

0.77 

(B-V) 

vs. 

(B-R) 

290 

0.78 

(B-V) 

vs. 

(R-D 

246 

0.44 

(B-V) 

vs. 

(V-D 

256 

0.61 

(B-V) 

vs. 

(V-R) 

292 

0.54 

(B-V) 

vs. 

Vc 

292 

-0.37 

(B-V) 

vs. 

i 

292 

-0.35 

Grt 

vs. 

(B-D 

228 

0.78 

Grt 

vs. 

(B-R) 

292 

0.76 

Grt 

vs. 

(B-V) 

292 

0.61 

Grt 

vs. 

(R-D 

279 

0.52 


+u. 


- 0 . 


.07 


).08 


-0.08 


+ 0 . 


continued on next page 


-0. 


FDR correction^ 


P-value^ 

CL(cry 

D 

i 2.5 cr' 

m 

^3crj 

m 

P 132.5c/ 

PI33ct‘ 

< 10-^ 

> 8.0 

25 

0.003690 

0.000804 

0.26 

0.30 

< 10-® 

> 8.0 

26 

0.003838 

0.000836 

0.25 

0.29 

< 10-® 

> 8.0 

27 

0.003986 

0.000868 

0.24 

0.27 

< 10-® 

> 8.0 

28 

0.004133 

0.000900 

0.26 

0.29 

< 10-® 

> 8.0 

29 

0.004281 

0.000932 

0.25 

0.28 

< 10-® 

> 8.0 

30 

0.004429 

0.000964 

0.23 

0.25 

< 10-® 

> 8.0 

31 

0.004576 

0.000996 

0.23 

0.25 

< 10-® 

5.42 

32 

0.004724 

0.001029 

0.23 

0.25 

< 10-® 

> 8.0 

33 

0.004871 

0.001061 

0.26 

0.29 

< 10-® 

> 8.0 

34 

0.005019 

0.001093 

0.23 

0.26 

< 10-® 

> 8.0 

35 

0.005167 

0.001125 

0.22 

0.25 

< 10-® 

> 8.0 

36 

0.005314 

0.001157 

0.23 

0.26 

< 10-® 

> 8.0 

37 

0.005462 

0.001189 

0.22 

0.25 

< 10-® 

> 8.0 

38 

0.005610 

0.001221 

0.22 

0.25 

0.000001 

4.89 

39 

0.005757 

0.001254 

0.26 

0.29 

0.000001 

4.83 

40 

0.005905 

0.001286 

0.24 

0.27 

0.000035 

4.14 

41 

0.006052 

0.001318 

0.23 

0.26 

0.000038 

4.12 

42 

0.006200 

0.001350 

0.23 

0.26 

0.000170 

3.76 

43 

0.006348 

0.001382 

0.22 

0.25 

0.000259 

3.65 

44 

0.006495 

0.001414 

0.22 

0.25 

0.000313 

3.60 

45 

0.006643 

0.001446 

0.23 

0.26 

0.000447 

3.51 

46 

0.006790 

0.001479 

0.26 

0.29 

0.000466 

3.50 

47 

0.006938 

0.001511 

0.23 

0.26 

0.000470 

3.50 

48 

0.007086 

0.001543 

0.26 

0.29 

0.000473 

3.50 

49 

0.007233 

0.001575 

0.23 

0.25 

0.000588 

3.44 

50 

0.007381 

0.001607 

0.23 

0.25 

0.001711 

3.14 

51 

0.007529 

0.001639 

0.24 

0.27 

0.002341 

3.04 

52 

0.007676 

0.001671 

0.24 

0.27 

0.002384 

3.04 

53 

0.007824 

0.001704 

0.22 

0.25 

0.004010 

2.88 

54 

0.007971 

0.001736 

0.22 

0.25 

0.004778 

2.82 

55 

0.008119 

0.001768 

0.22 

0.25 

0.005972 

2.75 

56 

0.008267 

0.001800 

0.22 

0.25 

0.007545 

2.67 

57 

0.008414 

0.001832 

0.24 

0.27 


< 10-® 

> 8.0 

1 

0.000148 

0.000032 

0.25 

0.28 

< 10-® 

> 8.0 

2 

0.000295 

0.000064 

0.25 

0.28 

< 10-® 

5.06 

3 

0.000443 

0.000096 

0.25 

0.28 

< 10-® 

5.22 

4 

0.000590 

0.000129 

0.22 

0.25 

< 10-® 

> 8.0 

5 

0.000738 

0.000161 

0.22 

0.25 

< 10-® 

> 8.0 

6 

0.000886 

0.000193 

0.22 

0.25 

< 10-® 

> 8.0 

7 

0.001033 

0.000225 

0.25 

0.28 

< 10-® 

> 8.0 

8 

0.001181 

0.000257 

0.22 

0.25 

< 10^® 

> 8.0 

9 

0.001329 

0.000289 

0.24 

0.27 

< 10-® 

> 8.0 

10 

0.001476 

0.000321 

0.24 

0.27 

< 10-® 

> 8.0 

11 

0.001624 

0.000354 

0.22 

0.25 

< 10-® 

> 8.0 

12 

0.001771 

0.000386 

0.22 

0.25 

< 10-® 

> 8.0 

13 

0.001919 

0.000418 

0.22 

0.25 

< 10-® 

> 8.0 

14 

0.002067 

0.000450 

0.25 

0.28 

< 10-® 

> 8.0 

15 

0.002214 

0.000482 

0.22 

0.25 

< 10-® 

> 8.0 

16 

0.002362 

0.000514 

0.22 

0.25 

< 10-® 

> 8.0 

17 

0.002510 

0.000546 

0.23 

0.26 
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Table 2. continued. 


FDR correction" 


x” 

VS. 

yc 

n"' 


Grt 

vs. 

Tn 

325 


Grt 

vs. 

(V-I) 

289 


Grt 

vs. 

(V-R) 

325 


Grt 

vs. 

Vc 

325 

- 

Grt 

vs. 

i 

325 

- 

Hr 

vs. 

a 

323 

- 

Hr 

vs. 

q 

323 

- 

iR-I) 

vs. 

(B-I) 

225 


{R-I) 

vs. 

(B-R) 

247 


iR-I) 

vs. 

(V-I) 

281 


(V-I) 

vs. 

(B-I) 

228 


iV-D 

vs. 

(B-R) 

257 


(V-I) 

vs. 

Vc 

309 

- 

(V-I) 

vs. 

i 

309 

- 

(V-R) 

vs. 

(B-I) 

228 


(V-R) 

vs. 

(B-R) 

292 


(V-R) 

vs. 

(R-I) 

280 


(V-R) 

vs. 

(V-D 

290 


(V-R) 

vs. 

Vc 

328 

- 

(V-R) 

vs. 

i 

328 

- 

(B-V) 

vs. 

Tn 

292 


(V-I) 

vs. 

Tn 

309 


(B-R) 

vs. 

q 

297 


Hr 

vs. 

e 

323 


(V-R) 

vs. 

Tn 

328 


Hr 

vs. 

Q 

323 

- 

(B-R) 

vs. 

e 

297 

- 

(B-V) 

vs. 

q 

292 


(B-I) 

vs. 

Tn 

228 


Hr 

vs. 

Vc 

323 


(B-I) 

vs. 

q 

228 


(B-V) 

vs. 

e 

292 

- 

(V-D 

vs. 

q 

309 


Grt 

vs. 

e 

325 

- 

Grt 

vs. 

q 

325 


(V-I) 

vs. 

e 

309 

- 

(B-D 

vs. 

e 

228 

- 

(R-I) 

vs. 

Vc 

281 

- 

— 3 cr FDR correction cutoff — 


Hr 

vs. 

(B-R) 

281 

- 

(R-I) 

vs. 

i 

281 

- 

(V-R) 

vs. 

e 

328 

- 

(V-R) 

vs. 

q 

328 


(R-I) 

vs. 

q 

281 


— 2.5 cr FDR correction cutoff - 

— 

IXROs* exceptforCKBOs 



(B-I) 

vs. 

(B-R) 

154 


(B-V) 

vs. 

(B-I) 

156 


(B-V) 

vs. 

(B-R) 

183 


(B-V) 

vs. 

(R-D 

163 


(B-V) 

vs. 

(V-I) 

168 


(B-V) 

vs. 

(V-R) 

185 


Grt 

vs. 

(B-I) 

156 





0.28!"' 

0.74+ 

0.72+ 

0.33! 

0.34! 

0.36! 

0.51! 

0.58! 

0.47! 

0.56! 

0.75! 

0.70! 

0.31! 

0.30! 

0 . 66 ! 

0 . 68 ! 

0.34! 

0.61! 

0.29! 

0.31! 

0.28! 

0.26! 

0.26! 

0.25! 

0.24! 

0.24! 

0.24! 

0.24! 

0.26! 

0 . 21 ! 

0.25! 

0 . 21 ! 

0 . 20 ! 

0 . 20 ! 

0 . 20 ! 

0.19! 

0 . 21 ! 

0.19+ 

0.17!°' 

0.17+ 

0.15+ 

0.15!; 

0.16+, 


.07 


continued on next page 


0 . 86 !“' 

0 . 79 !; 

0 . 81 !; 

0.46!; 

0.64!; 

0 . 61 !; 

0.81+, 


04 


07 


P-value^ CL{o-y 


< lo-** 

5.12 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

5.42 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10-® 

> 8.0 

< 10^® 

5.42 

< 10-® 

> 8.0 

0.000001 

4.85 

0.000005 

4.56 

0.000006 

4.52 

0.000006 

4.54 

0.000008 

4.47 

0.000016 

4.31 

0.000038 

4.12 

0.000040 

4.11 

0.000056 

4.03 

0.000117 

3.85 

0.000163 

3.77 

0.000294 

3.62 

0.000329 

3.59 

0.000350 

3.58 

0.000399 

3.54 

0.000667 

3.40 

0.001104 

3.26 

0.001697 

3.14 

0.004567 

2.84 

0.005277 

2.79 

0.005446 

2.78 

0.006765 

2.71 

0.007999 

2.65 


< 10 -® 

> 8.0 

< 10 -® 

> 8.0 

< 10 -® 

> 8.0 

< 10 -® 

> 8.0 

< 10 -® 

> 8.0 

< 10 -® 

> 8.0 

< 10 -® 

> 8.0 


D 

i 2.5 cr' 

m 

I-3crj 

m 

P02.5a-‘ 

18 

0.002657 

0.000579 

0.21 

19 

0.002805 

0.000611 

0.22 

20 

0.002952 

0.000643 

0.21 

21 

0.003100 

0.000675 

0.21 

22 

0.003248 

0.000707 

0.21 

23 

0.003395 

0.000739 

0.21 

24 

0.003543 

0.000771 

0.21 

25 

0.003690 

0.000804 

0.25 

26 

0.003838 

0.000836 

0.24 

27 

0.003986 

0.000868 

0.23 

28 

0.004133 

0.000900 

0.25 

29 

0.004281 

0.000932 

0.24 

30 

0.004429 

0.000964 

0.21 

31 

0.004576 

0.000996 

0.21 

32 

0.004724 

0.001029 

0.25 

33 

0.004871 

0.001061 

0.22 

34 

0.005019 

0.001093 

0.23 

35 

0.005167 

0.001125 

0.22 

36 

0.005314 

0.001157 

0.21 

37 

0.005462 

0.001189 

0.21 

38 

0.005610 

0.001221 

0.22 

39 

0.005757 

0.001254 

0.21 

40 

0.005905 

0.001286 

0.22 

41 

0.006052 

0.001318 

0.21 

42 

0.006200 

0.001350 

0.21 

43 

0.006348 

0.001382 

0.21 

44 

0.006495 

0.001414 

0.22 

45 

0.006643 

0.001446 

0.22 

46 

0.006790 

0.001479 

0.25 

47 

0.006938 

0.001511 

0.21 

48 

0.007086 

0.001543 

0.25 

49 

0.007233 

0.001575 

0.22 

50 

0.007381 

0.001607 

0.21 

51 

0.007529 

0.001639 

0.21 

52 

0.007676 

0.001671 

0.21 

53 

0.007824 

0.001704 

0.21 

54 

0.007971 

0.001736 

0.25 

55 

0.008119 

0.001768 

0.23 

56 

0.008267 

0.001800 

0.23 

57 

0.008414 

0.001832 

0.23 

58 

0.008562 

0.001864 

0.21 

59 

0.008710 

0.001896 

0.21 

60 

0.008857 

0.001929 

0.23 


1 0.000148 

2 0.000295 

3 0.000443 

4 0.000590 

5 0.000738 

6 0.000886 

7 0.001033 


0.000032 

0.30 

0.000064 

0.30 

0.000096 

0.28 

0.000129 

0.29 

0.000161 

0.29 

0.000193 

0.28 

0.000225 

0.30 


PyS3o-* 

0.24 

0.25 

0.24 

0.24 

0.24 

0.24 

0.24 

0.28 

0.27 

0.25 

0.28 

0.27 

0.24 

0.24 

0.28 

0.25 

0.25 

0.25 

0.24 

0.24 

0.25 

0.24 

0.25 

0.24 

0.24 

0.24 

0.25 

0.25 

0.28 

0.24 

0.28 

0.25 

0.24 

0.24 

0.24 

0.24 

0.28 

0.25 

0.25 

0.25 

0.24 

0.24 

0.25 


0.34 

0.34 

0.31 

0.33 

0.33 

0.31 

0.34 
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Table 2. continued. 


FDR correction" 


X’’ 

VS. 

yc 

n"' 

Grt 

vs. 

(B-R) 

183 

Grt 

vs. 

(B-V) 

185 

Grt 

vs. 

(R-I) 

175 

Grt 

vs. 

(V-I) 

180 

Grt 

vs. 

(V-R) 

197 

Hr 

vs. 

q 

191 

{R-I) 

vs. 

(B-I) 

154 

{R-I) 

vs. 

(B-R) 

162 

{R-I) 

vs. 

(V-I) 

176 

(V-I) 

vs. 

(B-I) 

156 

(V-I) 

vs. 

(B-R) 

166 

(V-R) 

vs. 

(B-I) 

156 

(V-R) 

vs. 

(B-R) 

183 

(V-R) 

vs. 

(V-I) 

180 

(V-R) 

vs. 

(R-I) 

175 

(B-R) 

vs. 

Vc 

186 

(B-R) 

vs. 

i 

186 

(B-R) 

vs. 

Tn 

186 

Grt 

vs. 

i 

197 

(B-V) 

vs. 

i 

185 

(B-V) 

vs. 

Vc 

185 

(V-R) 

vs. 

i 

197 

(B-V) 

vs. 

Tn 

185 

Grt 

vs. 

Vc 

197 

Grt 

vs. 

Tn 

197 


— 3 cr FDR correction cutoff — 


(B-I) 

VS. 

Tn 

156 

(B-I) 

vs. 

Vc 

156 

(V-R) 

vs. 

Vc 

197 

(V-I) 

vs. 

i 

182 

(V-I) 

vs. 

Vc 

182 

(B-I) 

vs. 

i 

156 

(V-I) 

vs. 

Tn 

182 

(V-R) 

vs. 

Tn 

197 

Hr 

vs. 

a 

191 

Hr 

vs. 

(B-R) 

176 

— 2.5 cr FDR correction cutoff 



P-value^ 

CL(cry 


i 2.5 cr' 

m 

^3crj 

m 

P 132.5c/ 

P/330- 

0.79-"° 

0.68!“ 

0.53+° 

05 

8? 

8? 

< 10° 

< 10-° 

< 10-° 

> 8.0 
> 8.0 
> 8.0 

8 

9 

10 

0.001181 

0.001329 

0.001476 

0.000257 

0.000289 

0.000321 

0.28 

0.28 

0.28 

0.31 

0.31 

0.32 

0.76+° 

8§ 

< 10-° 

> 8.0 

11 

0.001624 

0.000354 

0.28 

0.32 

0.75!° 

-0.39+° 

8§ 

06 

< 10-° 

> 8.0 

12 

0.001771 

0.000386 

0.27 

0.30 

07 

< 10-° 

> 8.0 

13 

0.001919 

0.000418 

0.27 

0.31 

0.62+° 

8^ 

< 10-° 

> 8.0 

14 

0.002067 

0.000450 

0.30 

0.34 

0.49+° 

8^ 

< 10-° 

> 8.0 

15 

0.002214 

0.000482 

0.30 

0.33 

0.59+° 

8? 

< 10-° 

> 8.0 

16 

0.002362 

0.000514 

0.28 

0.32 

0.78+° 


< 10-° 

> 8.0 

17 

0.002510 

0.000546 

0.30 

0.34 

0.72+° 

8? 

< 10-° 

> 8.0 

18 

0.002657 

0.000579 

0.29 

0.33 

0.69+° 

8§ 

< 10-° 

> 8.0 

19 

0.002805 

0.000611 

0.30 

0.34 

0.72+° 

8 ? 

< 10-° 

> 8.0 

20 

0.002952 

0.000643 

0.28 

0.31 

0.63+° 

8 ? 

< 10-° 

> 8.0 

21 

0.003100 

0.000675 

0.28 

0.32 

0.36+° 


0.000001 

4.85 

22 

0.003248 

0.000707 

0.28 

0.32 

-0.30+° 


0.000030 

4.17 

23 

0.003395 

0.000739 

0.28 

0.31 

-0.30+° 

8 ? 

0.000042 

4.09 

24 

0.003543 

0.000771 

0.28 

0.31 

0.29+° 

8 ? 

0.000076 

3.96 

25 

0.003690 

0.000804 

0.28 

0.31 

-0.28+° 

8 ? 

0.000076 

3.96 

26 

0.003838 

0.000836 

0.27 

0.30 

-0.27+° 

8 ? 

0.000194 

3.73 

27 

0.003986 

0.000868 

0.28 

0.31 

-0.27+° 

8 i 

0.000200 

3.72 

28 

0.004133 

0.000900 

0.28 

0.31 

-0.26+° 

8 i 

0.000262 

3.65 

29 

0.004281 

0.000932 

0.27 

0.30 

0.26+° 

8 i 

0.000271 

3.64 

30 

0.004429 

0.000964 

0.28 

0.31 

-0.25+° 


0.000468 

3.50 

31 

0.004576 

0.000996 

0.27 

0.30 

0.24!° 

8 ? 

07 

0.000792 

3.36 

32 

0.004724 

0.001029 

0.27 

0.30 

0.27+° 

08 

0.000796 

3.35 

33 

0.004871 

0.001061 

0.30 

0.34 

-0.27+° 

8 i 

0.000804 

3.35 

34 

0.005019 

0.001093 

0.30 

0.34 

-0.23+° 

8 i 

0.001228 

3.23 

35 

0.005167 

0.001125 

0.27 

0.30 

-0.24+° 

8 i 

0.001251 

3.23 

36 

0.005314 

0.001157 

0.28 

0.32 

-0.24+° 

8 i 

0.001263 

3.22 

37 

0.005462 

0.001189 

0.28 

0.32 

-0.25+° 

8 i 

0.001330 

3.21 

38 

0.005610 

0.001221 

0.30 

0.34 

0.22+° 

8 i 

0.002389 

3.04 

39 

0.005757 

0.001254 

0.28 

0.32 

0.21+° 

8 i 

0.002943 

2.97 

40 

0.005905 

0.001286 

0.27 

0.30 

-0.19+° 

8 ? 

0.007272 

2.68 

41 

0.006052 

0.001318 

0.27 

0.31 

-0.16!» 

8 ? 

08 

0.028693 

2.19 

42 

0.006200 

0.001350 

0.28 

0.32 


Computational steps of the correction for the false-discovery rate (see Sect. 3.41 X variable. Y variable. Number of data points. 
Spearman-p correlation value + the 68.2% uncertainty. ^ P-value of the correlation. 

Confidence level of the correlation in cr values of a normalized Gaussian distribution. 

Rank of the correlation after ordering by increasing P-values 
FDR-corrected P-value as a function of the number of tests for 2.5 cr conhdence. 

FDR-corrected P-value as a function of the number of tests for 3 cr conhdence. 

® Maximum detectable p at a 2.5 cr conhdence with a (3 risk of failing to detect it of 10%. 

Maximum detectable p at a 3 cr conhdence with a (3 risk of failing to detect it of 10%. 
















Table 5. Compilation of absolute magnitude and visible colors and corresponding orbital parameters 
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Object Clas^^ R^I VW B^ q Q i e a 7 ff References 

(42355) Typhon SDO(b) 7.36 + 0.08 13.97 + 3.01 0.78 + 0.11 0.52 + 0.04 0.42 + 0.12 0.93 + 0.13 1.71 +0.13 1.29 + 0.07 17.514 58.350 2.4 0.5383 37.932 2.61 0.0048212 0.540 2.684 4.684 19, 20, 35, 
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